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Abstract 


Gamma ray line emission from nuclear deexcitation following ener- 
getic particle reactions is evaluated. The con^iled nuclear data and 
the calculated gamma ray spectra and Intensities can be used for the 
study of astrophysical sites which contain large fluxes of energetic 
protons and nuclei. A detailed evaluation of gamma ray line production 
in the interstellar meditnn is made in the present paper. 
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I. INTRODUCTION 

The Interaction of energetic particles with ambient • matter produces 
nuclear gamma rays of energies ranging from tens of keV to about 20 MeV. 

These gamma rays exhibit a great wealth of spectral structure, ranging 
from, very narrow lines to broad features, depending on the composition and 
energy spectrum of the energetic particles, and the composition and physical 
state of the ambient medium. 

Observable gamma ray- line emission is expected from many astrophysical 

9 

sites, including solar flares, the interstellar medium, neutron stars and 
black holes, supernova remnants and nuclei of galaxies. Observations of such 
lines cai( reveal' directly many astrophysical processes so far studied only 
indirectly. Studies of diffuse gamma-ray line emission may give information 
on the intensity, spatial distribution and perhaps even the sources of the 
otherwise unobservable low-energy component (less than about 100 MeV/nucleon) 
of the cosmic rays. Gamma-ray spectroscopy may also help detemine the compo- 
sition of the interstellar medium on a galactic scale, independent of the ion- 
ization and molecular state of the material. Furthermore, very narrow lines 
with a width of a few keV at energies .of a few MeV permit the first direct 
measurements of the composition, size and spatial distribution of interstellar 
dust grains. 

Nuclear gamma rays were, in fact, observed from the solar flares of 1972, 
August 4 and 7 (Chupp. et al. 1973, 1975), and 1977, November 22 (Chambon 
et al. 1978), from the galactic center (Haymes et al. 1975, Leventhal, MacCallum 
and Stang, 1978a, b), from Centaurus A (Hall et al. 1976), and from a tran- 
sient event of about 20-minute duration observed from the general direction 
of the galactic anticenter (Jacobson \et al. 1978). 

Energetic particle interactions with ambient matter can lead to gamma 
ray line emission in a variety of ways. These include the direct excitation 
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of nuclear levels, the production of excited secondary nuclei and radioactive 
species which decay into excited levels of other nuclei, and the production 
of neutrons and positrons. Many excited levels decay by photon emission, 
and these gamma- rays are the topic of the present paper. 

Gamma ray line emission from radiative capture of neutrons and annihila- 
tion of positrons produced in energetic particle interactions has already 
been discussed in detail (Ramaty, Kozlovsky and Lingenfelter 1975) , and the 
positron annihilation line will be considered here for purposes of con^tarison 
only. 

In the present paper we provide a detailed treatment of gamma ray 
production from nuclear deexcitation following energetic particle inter- 
actions. To this end we have surveyed in detail the most important nuclear 
reactions of protons and alpha particles with the abundant constituents 
of cosmic matter (He, C, N, 0, Ne, Mg, A-t, - Si, S, Ca and Fe) , and we give 
the cross sections relevant for gamma ray line astronomy in Section II. 

The shapes of gamma ray lines contain valuable information on the 
physics of the emitting region, and they are important for considerations 
regarding the detectability of the lines above background. We treat gamma 
ray line shapes in Section III. Lines produced in the interstellar gas 
are broadened by the recoil velocity of the excited nucleus, but lines 
produced in the interstellar grains can be very narrow (Lingenfelter and 
Ramaty 1977) because some of the excited nuclei stop in solid materials 
before emitting gamma rays . We present several calculated gamma ray spectra, 
obtained by using a Monte Carlo simulation and the cross sections and line 
broadening effects . We also give the emissivitles of the strong lines 
normalized to the energy density and the energy deposition rates of the 
energetic particles, and we compare the production rate of positrons with 
those of the strong nuclear lines. The positron production rates are based 



on updated values of the cross sections given by Ramaty et al. (1975), 
which will be published Ip.^ a separate paper. 

The results of Sections II and III can be applied to the many astro- 
physical- sources that were mentioned above. Such -applications will become 
necessary with the expected launching in the near future of new gamma ray 
spectroscopy instrviments on satellites such as the Solar Maximum Mission 
(SMM) and the third High Energy Astrophysical Observatory (HEAO-C) . The 
general scope of gamma ray line astronomy, as well as a review of future 
missions, have been given recently (Lingenfelter and Ramaty 1978, Cline 
and Ramaty 1978). 

In Section IV we apply our results to gamma iray line emission from 
the interstellar medium. Previous treatments of gamma ray line production 
by low energy cosmic ray interactions with the interstellar gas were given 
by Rygg and Fishman (1973), Meneguzzi and Reeves (1975) and Lingenfelter 
and Ramaty (1976, 1977), using much less detailed nuclear data and a more 
limited choice of energetic particle spectra and compositions, and ambient 
medium compositions . In the present paper, we calculate the flimces of 
nuclear gamma ray lines from the interstellar .medium by taking into accotint 
all the available nuclear data, and we consider a broad range of spectral 
parameters and compositions. In particular, we evaluate the effects on 
nuclear gamma ray production of the possible spatial gradients of the relative 
abundances of heavy elements in galaxies (Peimbert, Torres-Peimbert and • 

Ray 1978, and references therein). We also compare the results of our 
calculations with observations of nuclear lines from the general direction 
of the galactic center. 


We summarize our results in Section V. 



II. NUCLEAR LINES AND CROSS SECTIONS 


In this section we presait the cross sections for the principal nuclear 
gamma ray lines resulting from interactions of energetic particles with 
ambient matter, assuming that both the particles and the matter have 
roughly solar composition. These lines include those produced both by 
direct excitation- and by spallation reactions leading either to nuclei in 
excited states or to radionuclei which decay to excited states. The esti- 
mation of which lines may be important depends on the cross sections for 
excitation of a level and the assumed elemental and isotopic abundances. 

The relative line intensities, which also depend of course on the energy 
spectrum of the energetic particles, will be calculated in section III. 

The lines which we have considered are listed by energy in Table I 
together with .the emission mechanism, the principal production processes 
and the meanlife of either the excited state or’ the radioactive parent if it 
is longer. The meanlife can be important in determining the line width as 
we shall discuss in section Illb below. For convenience the’ principal 
nuclei are also listed in Table 2 together with the relative solar abundance 
(Ross and Aller 1976), the direct deexcitation lines and the principal 
spallation product deexcitation lines . This' is the order in which the 
gamma ray line cross sections will be discussed. 
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a» Lines from He .Reactions 

Nuclear interactions of a-particles wxrn tie proauue t»ev«ici.j. 

important gamma-ray lines. The strongest of these are at 0.431 

' 7 7 

MeV and 0.478. MeV resulting from deexcitations in Be and Li, 
respectively (Tables 1 and 2) . The excited states of these isotopes 
C^Li* and ^Be*) can be populated by the reactions "^He (cc,n) ^Be* ,^He 
(ctfP)^l*i* and ^He Ca,n) ^Be (e) ^Li* (10%) , where the latter reaction 
includes production of Be in both the ground state and its excited 

. 7 

state. Since the mean life of Be is 76 days, decaying by electron 

7 

capture 10% of the time into Li*, the 0.478-MeV line resulting from 
7 7 - . 

Be(e) Li* is delayed in comparison with prompt lines such as those 
resulting from the direct excitation of ^Li* and ^Be*. ^Li and ^Be 
can also be produced from C, N and 0, but for compositions similar 
to solar abundances and steep energetic particle spectra these produc- 
tion -modes' ’are- quite negligible in comparison with aa reactions. 

The cross sections for aa reactions leading to total (ground' 
state plus first excited state) Be production and to Li production 
in its first- excited state are shown in Figure 1. The cross sections 

for the reactions ^He(a,p)^Li (total) and ^He (a,n) ^Be (total) were 

/ 

measured by King et aQ). (1975) from 9.75 Me.V/nucleon to 35 Mev/nucleon 
The cross section for the reaction "^He (a,p) ^Li* was measured from 
9.75 to 12.5 MeV/nucleon by S. M, Austin (private communication, 1975) 
Above 12.5 MeV/nucleon we assumed that the ratio of this cross section 

7 

to the total Li production cross section is the same as at 12.5 
MeV/nucleon. There are no direct measurements of the cross section 
for the reaction ^He (a,n) ^Be* ; we assume that above its threshold 
of 9... 5 MeV/nucleon, the cross section for this reaction is the 
same as that for the mirror reaction ^He (a>p) "^Li* . it 

7 

should be noted that the cross sections for Li (total) and 

7 

Be (total) are essentially equal (King et al. 1975) . A more 
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detailed discussion of these cross sections and their coinparison 

with previous estimates was given by Kozlovsky and Ramaty (1977) . 

The interactions of a-particles with He can also produce 

gamma xays at 35.562- MeV from ®Li*. As discussed by Kozlovsky 

and Ramaty (19 74b), because of isospin selection rules, this 

4 6 

level can be populated only by the reaction He(a,pn) Li*, and 

4 ^ 

not by He (a,d) Li*. For the former reaction we have adopted 

the cross section given by Mi tier (1972) , and we have assumed 

equal contributions for the production of 'Li* and °Li (g.s.). 

According to these estimates, the cross section for the reaction 
4 6 

He (a, pn) Li* is assiimed^ to be constant at about 1 mb above the 
threshold of about 15 MeV/nucleon. 

3 

Since in some astrophysical sites He may have large 

abundances (e.g. solar flares, Garrard, Stone .and Vogt. 1973), 

4 3 6 ‘ 

we also considered the reaction He (He, p) Li*. But since 
the cross section for this reaction, as measured by Harrison 
(1967) , is not large (< 10 mb) we do not take this processes into 
account' in our subsequent calculations. 

As discussed by Kozlovsky and Ramaty a974c), ^He has no excited 
states that decay primarily by gamma ray emission. However, 
photon deexcitation of some excited states is possible, but only 
states with isotopic spin T = 1 are of interest. The reason for 
this is that only such states can deexcite by dipole radiation, 
and hence have a reasonable chance of competing with particle 
emission. The energies of these states lie above 25 MeV (Fiarman 
and Meyerhof 1973) , bu-t exact values extracted from data on the 
reactions ^He(Y/P)^H, ^He(y,n)^He and ^H(p,n)^He and their in- 



7 


verses are quite model dependent (Gibson 1972) . Two possible 

levels may lie at ~27 MeV and ~30 MeV. Kozlovsky and-Ramaty 

\ 

(1974c) have estimated that the upper limits on the cross 

section for gamma-ray emission from these levels following (PfP') 

“3 

excitation is only about 10 mb. This value is quite low, and 

4 

hence gamma-ray lines from (prp') reactions on He do not appear 
to have important astrophysical consequences. 

b. Lines from C, N and O' Reactions 

12 

The strongest deexcitation line in C is at 4.438 MeV re- 
sulting from. the deexcitation of its first excited, state at 4.439, 

MeV. The various excitation modes of this level are given in 

12 12 4 439 

Table 1. The cross section for the reaction C(p,p') C* 
has been measured by Reich et al. (195 6)' between 5 and 5.7 

MeV, by, Barnard et al. (1966) between 6 and 11.5 MeV, by Conzett 

(1957) between 10 and 12 MeV, by Daehnick and Sherr (1964) between 

14 and 19 MeV, by Dickens et al. (1963) between 18 and 30 MeV, 

by Stovall and Hintz (1964) at 40 MeV, by Fannon et al. (1967) 

at 49.5 MeV, by Horowitz and Bell (1970) at 100 MeV, by Emmerson 

et al. (1966) at 145 MeV, and by Tyren^ and Maris (1957) at 185 

MeV. This cross section is. shown by the solid curve in Figure 2. 

We have used the average cross section in the resonances near 

the maximum at 10 MeV. However, we have not averaged the cross 

section of the resonance at 5.35 MeV since it may have some effect 

on the calculations for very steep particle spectra. 

The data points represented by the closed circles and squares 

are based. on measurements of gamma rays at 4.44 MeV resulting from 
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the bombardment of C with protons (Zobel et al. 1968, Alard et al. 
1974) . The closed circle points were deduced by using the differ- 
ential cross section at 135*^ given by Zobel et al. (1968) and 
their data oh the angular distributions of gamma rays produced 
in proton bombardment of C. The data point of Alard et al. (1974) 
is already given in integral form. The difference between this 
gamma ray data (dashed curve) and the (p,p') data represented by 
the solid line is 'most likely due to the excitation of the 4.444 
MeV level of B by the reaction C(p,2p) B* which has an energy 
threshold of 22 MeV. Even though the 4.438 and 4.443 MeV lines 
cannot be resolved because kinematical Doppler broadening blends 
them into a single feature, in our calculations we have separated 
these lines by taking the cross section for producing the 4.443 
MeV line from C equal to the difference between the dashed and 
solid curves in Figure 2. 

The cross sections for reaction ^^C(a,a' ) are 

given by the dashed-dotted curve in this figure. This curve is 
based on measurements of Mitchel, Carter and Davis (1964). between 
1.5 and 4.3 MeV/nucleon, of Corelli, Bleuler and Tendam (1959) 

.at 4.6 MeV/nucleon, of Yavin and Farwell (1959) at 10 MeV/nucleon. 
The open circle at 13 MeV/nucleon (Zobel et al. 1968) is data 
based on gamma ray measurements; thus, as discussed above, it 
includes' the contribution .of the reaction ^^C(a,a* ) 

The dashed-crossed curve is our estimate for the sum of the cross 
sections of the reactions ^^C (a,a' ) and ^^C(a,x) 

As in the case of proton induced interactions, wfe interpret the 
difference between the gamma ray data and the (a,a') data as due 
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to the contribution of In our calculations, we use the 

dashed-dotted curve up to 5 MeV/nucleon and'’ the dashed-crossed 
curve at higher energies. ‘As above, we separate the 4.438 and 
4.443 Me'V lines by taking the difference between the dashed-dotted 
and dashed-crossed curves. 

12 

The 4.439 MeV level in C can also be populated by spallation 

reactions on and .as indicated in Table 1. The cross 

16 

section's for reactions on 0 were taken from Zobel et al. (1968) 
and Alard et al. (1974) , and are shown by curve 1 in Figure 3 for • 
proton induced reactions, and by curve 6 for a particle induced 
reactions. The shape of curve 1 below 25 Me'V was obtained by 
normalizing the- measurements of Zobel et al. (1968) to the pre- 
liminary data of Dyer and Bodansky (private communication 1977) 

for the reaction ^'^Mg(p,3CYj^ 634^^^^®* Since both this and the 
16 12 

0(p,x) C* reactionary (p,.pa) processes, we expect them to 

" 16 12 

have similar energy dependences. For the reaction d'(a,xy^ 44 ) ^ 

we have only one measurement (Zobel et al. 1968) of -110 mb at 

13 Me'V/nucleon , and curve 6 is our estimate of the energy dependence 

of the cross section for this reaction. For proton spallation of 

14 

N we use the measurement of Clegg et al. (1961) at 120 MeV, from 

which we estimate that the 4.438 MeV gamma ray cross section from 

is larger than that from ^^0 by about a factor of 2.4. -With 

14 

solar abundances, N/0 ~0.13, the 4.438 MeV yield from N is 'about 

30% of the yield from 

. 12 ' 

Exerted states of C above the 4.439 MeV level decay mostly 
by particle emission and hence they are not important sources of 
gamma ray lines. An exception is. the 15.11 MeV level, which, be- 
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cause of conservation of isotopic spin, cannot decay fay emitting 
a particles and hence deexcites only by gamma ray emission. The 
cross sections for exciting this level were compiled by Crarinell, 
Ramaty and Crannell. (19-77)- who found that the intensity of the 

15.1 MeV line is at most about 2% of the 4.438 MeV line intensity. 

12 ^ 
Spallation reactions in C can produce several gajnma ray 

11 11 

lines (Table 2) . The first exerted states of C and B at 
1.995 MeV and 2.124 MeV can be populated by (p,pn) and (p,2p) 
reactions. Zobel et al. (1968) have measured the combined cross 
section for gamma-ray emission from these two levels and their 
data is shown in Figure 3 by open triangles. Curve 4 in this 
figure gives our estimate of this cross section, where we have 

, X 

used the fact that the thresholds of the above reactions are 
about 20 MeV. In our subsequent calculations, we take the cross 
sections of these two lines equal to each other. Zobel et -al.- 
(.1968) have also observed a cluster of four lines between ~6.3. 

and 6.8 MeV at the energies indicated in Table 1, due to the 

• . 11 11 
deexertatron of higher levels in C and . B. The cross section 

for this cluster, as given by Zobel et al. (1968) and Clegg et al. 
(1961) , is on the average about 25% of the cross section (Figure 3) 
for the production of the ~5.2-MeV feature from 0 by the same 
type of reactions. We use this ratio in our subsequent calcula- 
tions and we assume equal cross sections for the four lines, 

' 12 

Two other strong spallation lines from C are at 0.717 and 

10 * > 

1.023 MeV resulting from deexcitation in B. The cross sections 
for the production of these lines, measured by Clegg et al. (1961) 
at proton energies around 140 MeV, are, respectively, about 50% 
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and 20% of the cross section for producing the 4.44 MeV line from 
X6 

0 spallation at the same energy. Curve 3 in Figure 3 is our 

estimate of the cross section for 0.717 MeV photon production 
12 

from Cr where we use the fact that the threshold for this process 

is about 20 MeV. For the 1.023 MeV photon production cross section 

we take 40% of the values represented by curve 3 in Figure 3. 

13 

The principal lines from deexcitation of C are at 3.684’ and 

3.853 MeV {Tables 1 and 2). The cross section for proton excitation 

of these levels, measured by Guratzch et- al. (1969) at 7 MeV, are, 

respectively, about 80% and 30% of that for the 4.439 MeV level 
12 

xn ^ C. Using the branching ratios of de Meijer, Plendl and Hol\ib 

(1974), we estimate the 3.684 and 3.853 MeV gamma ray production 

cross sections from (p,p') reaction on to be, respectively, 

90% and 20% of that for 4.438 MeV gamma rays from (p,p') reactions 
12 

on C shown in Figure 2. For the production cross sections of 

these lines by (a, a') reactions on we assume that these 

ratios also hold with respect to the 4.438 MeV production cross 

12 

sectxons from (a,a’) reactions on C. 

The strongest line from the deexcitation of is at 2.313 

MeV. The various processes leading to photons at this energy 


are listed in Table 1, and their cr-dss sections are shown in 

Figure 4. The cross section for the reaction ^^N(p,p '7 )-^^N, 

2 . 31 3 

shown by the solid curve, includes- the deexcitations of all the 
. 14 « 

levels xn N whxch cascade to the ground state through the 2.313 


MeV level. The data from 3.8 to 6.4 MeV are from Phillips et al. 
(1972), from 6.5 to 10 MeV from Boreli et al. (1968), at 10.2 MeV 
from Donovan et al. (1964),- and from 9 to 26 MeV from Hansen et al. 
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(1973) . The branching ratios are from de Meijer et al. (1974) . 

The other two strong deexcitation lines from are at 1.632 

MeV and 5.105 MeV (Tables 1 and 2). The cross sections for the 

production of these lines.,, derived from the above measurements ■ 

and branching ■ ratios , are approximately 60% and 70% of the cross 

14 14 

section for the reaction N(p,p'Y 2 313 ^ ^ shown by the solid 

14 

curve in Figure 4 . Other lines resulting from cascades in N 

have cross sections less than 10% of the 2.313 MeV line cross section 

14 ' 14 

The -cross section for the reaction N(p,XY 2 ^ 33 ^ 3 ) N 
shown by dashed curve is the sum of the cross section shown 

by the solid curve and of the cross section of the reaction. ^"^N(p,n) ^^ 0 , 

. 14 ■ 1 4 *2 31 3 

since 0 decays essentially 100% of the time to N * . The 

cross ' section for this (p,n) reaction from 6.5 to 12 MeV is taken 

from Kuan and Risser (1964) . This cross section, however, is 

quite uncertain, since preliminary data from Dyer and Bodansky 

(private communication 1977) suggests that the cross section for 
the reaction ^^N(p,n)^^0 may be much lower than reported by Kuan 
and Risser (1964). The cross section for the reaction ^^N(p,xy 2 ^ 33 _ 3 j 
between 20 and 24 MeV is from the preliminary data of' Dyer and 
Bodansky (1977) and the data point at 120 MeV is from Clegg et al. 
(1961). 

The dashed-crossed curve in Figure 4 is the cross section 
for the reacion N(a,a*Y 2 3 x 3 ) “ cross section is based 

on the measurements of Ploughe (1961) at 4.8 MeV/nucleon, Garcia, 

Milio and Senent (1970) at 5.4 MeV nucleon, and Harvey et al. 

(1966) at 10.1 MeV/nucleon, and the branching ratios of de Meijer 
et al. (1974) . We assume that the relative intensities of the 
5.105 and the 1.632 MeV lines compared to the 2.313 MeV are the 
same for (a,a’) excitations as for^,p') excitations, i.e. about 
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70% and 60%, respectively. These values are consistent with the 

measurements of Garcia et al. (1970) . 

Gamma rays at 2.313 MeV resulting from proton interactions with ^^0 

(Tables 1 and 2) have been measured by Zobel et al. (1968) and their data is 

shown by the closed circles- in Figur^ 4. These gamma rays result from the 

^ 1 Ll 

direct production of N and from the decay of 0, with a possible contribu— 

13 

tion from excited states 'of N. The dashed-dotted curve through the data 

points of Zobel et al. (1968) is our estimate of the cross section for 2.313 

X6 

MeV photon production from 0. 

16 

Nuclear excitation of 0 can produce several gamma ray lines (Table 2) . 

16 + — 

The first excited state of 0 at 6.05 MeV decays by e - e pair emission. 

The second, third and fourth excited states at 6.131, 6.919 and 7.119 MeV 
decay almost exclusively- to the ground state producing gamma rays at 6.129, 
6.917 and 7.117 MeV, respectively. The fifth excited state at 8.872 MeV 
decays 75% of the time to 6.131 MeV level and hence produces lines at 2.741 
MeV and 6.129 MeV. Other transitions between the states of 0 are also 
possible (de Meijer et al. 1974). In particular, excited levels at 10.94 and 
11.07-MeV cascade to the ground state via the excited states at 6.131, 6.919 
and 7.119 MeV with branching ratios given by Ajzenberg-Selove and Lauritsen 
(1968). 

The cross section for producing gamma rays at 6.129 MeV from proton 
16 

interaction with 0 is shown in Figure 5 by the solid curve. This cross 

T ^ ^ ft 1 ”1 * 

section is sum of the cross sections of the reaction °0(p,p*) °0 ' and 

76% of the cross section for the reaction ^^0(p,p*) ^^0 and about 50% 

of the cross section of the reaction ^^0(p,p^) ^^0 These cross sections 

were measured by Dangle et al. (1964) between 7 and 10.5 MeV, by Kobayashi 
(1960) between 11 and 15 MeV, by Daehnick (1964) between 15 and 19 MeV, by 
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Homyak and Sherr (1955) at 19 MeV, by Crawley and Garvey (1967) at 17.5 MeV, 
by Austin et al. (1971) between 17 and 45 MeV, by Sundbeirg and Tibell (1969) 
at 185 MeV, and by Friedes et al. (1967) at 1 GeV. 

-The-dashed-dotted" curve in Figure 5 is the cross section for’ the reaction 

'’0(a,a’y^ i«n) ”0. This cross section was measured by Mehta, Hunt and Davis 
6.129 

(1967) from 2.5 to 4.7 MeV/nucleon, by Corelli, Bleuler and Tendam (1959) at 
4.5 MeV/nucleon, by Blatchley and Bent (1965) at 5.6 MeV/nucleon, by Yavin 
and Farwell(1959) at 10 MeV/nucleon, and by Hairvey et al. (1964) at 16,2 
MeV/nucleon. 

Most of the experimental measurements of the cross sections for excita- 

16 

tion of the 6.919 and- 7.119 MeV levels in 0 by (p,p') reactions do not 
resolve these two levels. Therefore, in Figure 6 the solid curve gives a fit 
to the sum of the cross sections for the excitation of these levels based on 
measurements of Zobel et al. (1968). Measurements of tp,p’) reactions to these 
levels by Kobayashi (1960) at 5.6 MeV and Crawley and Garvey (1967) at 17.5 - 
MeV coincide with this curve. At. higher energies the measurements of (p,p’) 
reactions by Hornyak and Sherr (1955) and Sundberg and Tibell (1969) give 
lower values. The difference should be attributed to contributions of spalla- 
tion lines such as the 7.299 MeV from which show up at high incident 

proton energies. In our calculations we use the cross section of Figure 6, 
and we asstime that half of the resultant photons are at 6.917 MeV and half are 
at 7.117 MeV, 

The 2.741 MeV line is produced by transitions between the 8.872 level 
to the 6.131 MeV level in ^^0 with a branching ratio 76%. The dashed curve 
in Figure 6 is based on the (p,p') measurements of Crawley and Garvey (19-67) 
at 17.3 MeV, and of Austin et al. (1971) between 17 and 45 MeV. The data of 
Zobel et al. (1968) are in agreement with the (p,p’) measurements.. The 
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^^ 0 (a, 3 cr^ Q - , -7 or ^^Ncxoss section is shown by the dashed- 

/ •«3 

/ 

dotted curve in Figure 6 . This curve is based on measurements' of 

Corelli et al. (1959) at 4.5 MeV/nucleon, of Harvey et al. (1966) 

at 10 MeV/nucleon and of Zobei et al. (1963) at 13 MeV/nucleon. 

1 

Gamma ray lines from 0 can also be produced by spallation 

12 15 15 14 

reactions which populate excited levels in C, N, O, O, . 

and (Table 2). The production of 4.438 MeV gamma rays 

16 12 

by the reaction 0(p,xY^ ^ 3 g) C has been discussed above. 

’* 15 

The strongest lines from N are at 6.322, 5.270 and 5.298 

MeV and from ^^0 at 5.180, 5.241 and 6.176 MeV. Because of poor 

energy resolution the 6.2 MeV gamma ray production cross sections of 

Zobelet al. (1968) and Alard et al. (1976) shown by the data points 

1 6 

in Figure 5 include the lines at 6.129 MeV from O, at 6.176 from 
O and 6.322 MeV from N. Thus, the difference between the dashed 
and solid curves, and between the dashed-crossed and dashed-dotted 
curves, in Figure 5 give the sum of the cross sections for the 
production of the 6.176 and 6.322 MeV lines in proton and alpha 
particle induced reactions, respectively. There is only one direct 
high-resolution measurement of gamma rays from ^^0 which can re- 
solve the 6.129, 6.176 and 6.322-MeV lines (Goryachev et al. 1973). 
This measurement at ~1 GeV shows the 6 . 322-MeV. line to be stronger 
than the line at 6.176 MeV by about a factor of 2. Assuming that 
this ratio holds also at lower energies, we take the cross sections 
for the production of the 6.322-MeV and the 6.176-MeV lines as 2/3 
and 1/3, respectively, of the combined cross sections of these two 


lines . 
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16 14 * 

The sum of the cross sections for the reactions 0(p,x) N , 

16 16 * 16 16 * ' 

0(p,^pn) 0 , and 0(p,2p) leading to ~5.2 MeV photons is 

shown by curve 2 in Figure 3 based on data by Zobel -et aL. .(.19.68) . 

In our calculations we arbitrarily assume that the cross sections 

for the production of the 5.105, 5.180, 5.241, 5.270 and 5.298 

MeV lines are each 20% of this cross section. The cross section 

for the production of ~5.2 MeV photons in a particle induced 

spallation reactions is shown by curve 7 in Figure 3. As with the 

proton induced reactions, we assume that the cross sections for the 

individual lines are 20% of the total. " ' 

The lines at 0.717 and 1.023 MeV can be produced by the 
16 

spallation of 0. Foley et al. (1962a) find that for 140 MeV 

protons the cross section of these two lines are about equal to 

each other, and ‘each of them is about half of the cross section 

for producing the 4.438 MeV line from 0 spallation. We therefore 

approximate each of these cross sections by curve 3 in Figure 3. 

Spallation of O can also populate excited levels of leading 

to iine emission at 3.684 MeV and 3.853 MeV. In Figure 3 we show 

> 

the data of Zobel et al. (1968) for the combined production of 
both these lines and our estimate (curve 5) for the energy depend- 
ence of this cross section. 

c. Lines from He,' Mg, Al # Si, S, and Fe Reactions 

n 

The strongest deexcxtation line from Ne is at- 1.634 MeV 
resulting from the deexcitation of its first excited state. The 
various production modes of this line are given in Table 1. The 
'solid curve ^in Figure 7 is the cross section for the reaction 
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20 20 *1 fi ’^4 

Ne(p,p') Ne * . From 5 to 14.2 MeV. the data are from Oda 

et al. '(I960) and references therein, at 17 MeV from Schrank et al. 
(1962), and at 24.5 MeV- from de Swiniarski et al. (1969). At low 

energy we have averaged the data over resonances. 

20 20 

The cross section for the reaction Ne(p,p'Yj^ 634^ which 

includes the cascades from the 4.247 and 4.968 MeV level's is shown 
by the dashed curve in Figure 7. The cross sections for populating 
these levels were measured by Oda et al. (1960) from 7,6 to 14.2 
MeV, Schrank et al. (1962) at 17 MeV, and de Swiniarski et al. 

(1969) at 24.5 MeV, and we have used the branching ratios of 
de Meijer at al. ’(1974). Since there is no data above. 24. 5 MeV, 
the dashed curve at these energies is an approximation based on 
similar excitation functions of ^^0 and ^^Mg. Below 5 MeV we have 
used preliminary data from gamma ray measurements (P. Dyer, and 
D. Bodansky, private communication 1977) which at higher energies 
aiB consistent with the above measurements and calculations. 

The dashed-dotted curve in Figure 7 is the cross section 

for the reaction ^Ne (a,a' ) * . This cross section has 

been measured by Seidlitz-, Bleuler and Tendam (1958) at 9.5 MeV/ 

nucleon and by Rebel et al. (1972) at 41 MeV/nucleon. There is 

no data on this cross section at lower energies and we have based 

our estimate on the expected behavior of the cross section near 

the Coulomb barrier. The dashed-crossed curve is the cross section 

20 20 

for. the reaction Ne(a,a*Y-L 534 ) N®* calculate 

the contributions of the cascades from the measurements of Seidlitz 
et al. (195'8) , and we assume that the* ratio ("2) between the cross 
sections for the reactions ^°Ne(a,a'Yj^ g 2 ^)^°Ne and ^®Ne (a,a' ) ^^Ne*^' 
remains constant at higher energies. 

Transitions from the 4.247 and 4.968 MeV levels provide the 



18 


20 

next strongest deexcitation lines in Ne, at 2,613 MeV and 3.334 

MeV. At 17 MeV Schrank et al. (1962) find that the intensities 

of these lines are about 25% and 10% of the 1.634 MeV line. We 

have assumed that these- ratios are the same at all energies. 

20 

There is no data on spallation gamma ray lines from Ne. 

But as can be seen from the measurements of Zobel et al. (1968) 

and Chang et al. (1974) the most important proton- induced reactions 

20 

are (p,n ) , (p,2p) , (p,pn) and P/pa) . For Ne (P/n) reactions 

20 . 20 
lead to Na which is radioactive and decays back to Ne producing 

gamma rays at 1.634 MeV about 80% of the time. Because of the short 

. 20 . s . 

mean life (0..64 sec) of Na, the contributions of these gamma rays 
are included in the measurements of Dyer and Bodansky (private 
cummunication- 1977) . Prom the comparison of these measurements 
and the (p,.p') data discussed above, we estimate that the contri- 
bution of Na to the total 1.634 MeV line is less than a few 
percent. For the other reactions, the strongest lines are likely 
-to result from the deexcitation of the first two excited states of 

the product nuclei. Therefore, we expect lines at 0.110 MeV and 

19 19 

0.197 MeV from P at 0.238 MeV and 0.275 MeV from Ne and at 

X6 

6.129 MeV from 0. We assume that the sum of the cross sections of 

19 19 

the four lines from F and Ne equals the cross section for 

~5.2 MeV photon production from spallation given by curve 2 

in Figure 3, and that the cross section for the reaction 

^^Ne(p,XYg 129^^^^ same as the cross section for the reaction 

16 12 

0(P/xy^ 438' ^ given by curve 1 in this figure. 

OO 

' The principal deexcitation line from Ne is at 1.275 NeV 
from the first level. Measurements of the cross section for proton 
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excitation of this level have been made at 1.9 to 3.2 MeV by Sorokin 

et al. (19 6'3) , at 4.8 to 14 MeV by Hulubei et al. (1969), and at 

24.5 MeV by de Swiniarski et al. (1972). These measurements suggest 

that this cross section is roughly 2/3 of that for proton excitation 

of the 1-634 MeV level in ^^Ne. We therefore assume that the ^^Ne 

(P,P*Yj^ 275^ Ne(a,a'Y, 2 ^ 27 cross sections are equal 

20 20 

to the product of this fraction and the Ne(p,p'Y]^ 634^ 

20 20 

the Ne(a,a'Y 2 534 ) cross sections, respectively, in Figure 7. 

The nuclear interactions between energetic particles and comple 
nuclei (A^20) produce many more gamma ray lines than those result- 
ing from the deexcitation of the low-lying levels discussed above. 
These lines result from the cascade deexcitations of the many high- 
lying level's that are populated both by direct excitations and 
spallation reactions. Zobel et al. (1968) have measured the total 
production cross section, 0 ^, of gamma rays of energies greater 
than 0.7 MeV in interactions of protons and alpha particles with 
complex nuclei. They find that for C and 0 this total cross section 
is essentially the same as the sum of the cross sections for the 
excitation of individual lines; but for heavier targets the resol vabl 
lines can account for only a fraction of the total gamma ray 

cross section, 0 . 

Y 

The measurements of Zobel et al. (1968) for 0 ^ were done, at 
proton energies of 16 MeV, 33 MeV, 56 MeV and 160 MeV for C, O, 

Mg, Al, Fe and Co targets and at an a particle' energy of 15 MeV/ 
nucleon for C, 0, Al and Fe targets . We have interpolated these 
results for Ne, and the resultant cross sections at the above 
energies are shown fay the diamonds in Figure 7. The crossed and 
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crossed-dotted curves are our estimates of tJie energy dependences 
of for proton and a particle induced reactions, respectively. 

At high energies we have assvimed a constant ratio between these 
two curves, and at low energies we have merged them with the ap- 
propriate cross sections for the excitation of the first excited 
level of ^^Ne. 

In addition to the lines discussed above, gamma ray emission 

is also produced from the decay of long lived radioactive spallation 

products. Radioisotopes with mean lives longer than 1 day that we consider are 

are ^Be, ^^Na, ^^Co and ^^Co. We have already 

• - 7 

discussed the production of Be in Ila, In order of increasing 
nuclear mass, we next consider Na which produces a gamma ray line 
at 1.275 NeV by decaying into the first excited- state of ^^Ne. 

, 09 

For Ne isotopes, an important production mode of Na is 

the reaction ^^Ne (p,n) ^^Na. Although the cross section of this 

reaction has not been measured, we estimate from the systeraatics 

.of other (p,n) reactions that this cross section should be roughly 

26 26 

equal to that of Mg{p,n) Al (shown in Figure 9) and we use 
these values in our calculations. The production of ^^Na from 

Mg* and Si is discussed below. 

1 

The strongest deexcitation line in ^'^Mg is at 1.369 MeV. 

The solid curve in Figure 8 is the cross section for the reaction 
24 24 *1 369 

Mg(p,p') Mg ' . -The data for this reaction from 3 to 5,5 

MeV is from Duray et al. (1972), at 5. 4- to 7 MeV from Seward (1959), 
at 12 MeV from Conzett (1957), at 17.5 MeV from Crawley and Garvey 
(1967) , at 40 MeV from Stoval and Hints (1964) , and at 100 MeV 
from Horowitz et al. (1969). 
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The cross section for gamma ray emission at 1.369 MeV from 
proton bombardment of including cascades from higher levels, 

has -been measured .at 30.3 MeV by Zobel et al. (1968) and at 14 3 
MeV by Foley et al. (1962b) . These data are shown by the open 
square and triangle in Figure 3. These measurements are consistent 
with estimates of the total 1.369 MeV emission mode from the meas- 
urements of Crawley and’ Garvey (1967) and Horowitz et al. (1969), - 
* 

using the branching ratios given by de Meijer, Drentje and Plendl 

(1975) and Endt and Van der Leun (1973)'. These estimates suggest 

that the contribution of the cascades from the first 10 excited 

states out of about 40 bound states is about equal to -the contribution 

of the direct excitation of the first excited levelJ The dashed 

curve in Figure 8 is our estimate of the cross section for the 
24 24 

reaction Mg(p,p'y^ 369*^' find that thxs estimate is in 

good agreement with the preliminary gamma ray data of Dyer and 
Bodanskyf (private communication 1977) . 

The cross section for the reaction ^^Mg(a,a' ) ^"^Mg is 

shown by the dashed-dotted curve in Figure 8. This cross ection 

i 

has been measured by Sauter arid Singh (19 74) from -2.5 to 3.5 MeV/ 
nucleon, by Eberhard and Trombik (1972) from 3.9 to 4.8 MeV/nucleon, 
by McNeilley et al. (1973) at 4.2 MeV/nucleon, and by Janus and 
McCarthy (1974) at 10 MeV/nucleon. Since there is no data below 
2.5 MeV/nucleon, our estimate in -this energy range is based on 
the expected behavior of the cross section near -fche Coulomb bar- 
rier. The -1.369 MeV level will also be populated 'by cascading 
from highler lying" levels . Prom Eberhard and Trombik (19 72), who 
have measured -the excitation cross sections of the 4.122, 4.239 
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26 

and 5.236 MeV leve.ls, we estimate that the cascades in Mg 
increase the effective 1.369 MeV gamma- ray production cross 
section by about a factor of 2.5. Similarly, from the measure- 
ments of Thompson et. al. (19-6 7) at 2.7 MeV/nucleoh , we find that 

this factor is about 1.4. The dashed-crossed curve in Figure 8 

24 

is our estimate of the cross section for the reaction Mg(a,a'*Y^ 369^ 

2 %. 

Although many more deexcitation lines can result from transi- 

24 

trons among the 40 bound states of Mg, we estimate from the 

cross sections of Crawley and Garvey (1967) and the branching ratios 

of.de Meijer et al. C1975) that none of these lines is stronger 

than 10% of the 1.369-MeV line. These lines are included in the 

24 

total unresolved nuclear radiation from Mg. In Table 1 we include 

only the Tine at 2.754 MeV. This line results from transitions 

between the second and first excitated. states of ^"^Mg, and estimate 

that its intensity is 10% of the T.369 MeV line intensity.' 

The two strongest lines resulting from the spallation of 
24 

Mg are at 1.634 and 1.636 MeV. These lines are due to the de- 
. 20 23 

excitation of Ne and Na which result from '(p,pa) and (p,2p) 

reactions, respectively. Zobel et al, (1968) have observed a 

feature at ~ 1.6 MeV from the bombardment of ^^Mg with 30 MeV 

protons, and found that its intensity is about twice that of the 

16 

4.44 MeV- line from 0, whereas Foley et al. (1962b) found that 

this ratio at 143 MeV is about unity. In our calculations we 

tcike the sum of the cross section for producing 1.6 34 and 1.636 

24 

MeV photons from Mg spallation to be,' at all incident energies, 

1 • 12 

1.5 times the. cross section for the reaction °0(p,XY^ 438-^ ^ 
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shown in Figure 3. This assumption is consistent with the prelim- 
inary data of Dyer and Bodansky (private communication 1977) for 
incident proton energies from about 16 to 23 MeV. We further 
assume that these lines have equal cross sections. The 1.6 MeV 

feature observed by Zobel et al. (1968) could also include the 

23 

1.600 MeV line of Mg, and- the comparison, between these measure- 
ments and those of Dyer and Bodansky (which have much better energy 

> 

resolution) suggests that the cross section for producing 1.600 
MeV photons from ^^Mg spallation is about 50% of the ^^0(p,xy^ 438^ 
cross section. 

Other spallation lines from ^^Mg are at 2.613 MeV and 2.640 
MeV resulting from deexcitatiohs in ^^Ne and ^^Na, respectively. 
From the measurements of Zobel et al. (1968) , at an incident proton 

energy of 30 MeV, the sum of the cross sections of these two lines 

16 12 ' * 
is about 60% of the 0(p,xy^ 438^ section, whereas 

from the measurements of Foley et al. (1962b) , at 143 MeV the 

corresponding ratio is about 1.5. '‘We assume a constant ratio of 2, 

and we take the cross sections for producing the 2.613* and 2.6 40 

MeV lines equal to ea!ch other. Foley et al. (1962) observed , 

another spallation feature at 0.92 MeV with half the strength 

of the 4.438 MeV line from at ~ 143 MeV. They identified 

22 

this feature as the 0.891 MeV line from Na, and we take its • 
cross section, at all incident energies, equal to 0.5 times the 
0(p,xy^ 438^^^*" cross section. 

Neither Foley et al. (1962b) nor Zobel et al, (1968) looked 
for gamma-rays below 0.7 MeV, but both the 1.636 and 1.600 MeV 
lines imply additional lines of at least equal intensities, at 
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0.440 and 0.451 MeV, respectively. 

24 

Because of the multitude of excited states of Mg and its 

spallation products, nuclear interactions produce many weaker 

lines in addition to the lines that we have discussed above. 

The crossed and crossed- dotted curves in Figure ‘8 show the total 

24 

gamma ray production in Mg. The closed square is a measurement 

by Zobel et al. (1968) , while the diamonds have been obtained' from 

interpolations similar to those done for Ne. At low incident 

energies, these total cross sections approach those for the 

direct excitation of the i.369 MeV level. 

22 

The cross section for' Na production from Mg is shown in 

Figure 9. Measurements of -tdiis cross section for proton spalla- 
24 • 25 

tion of bo-fh Mg and Mg are reported and reviewed by Furukawa 

et al. (1971) for proton energies up to 52 MeV, by Korteling 

and Caretto (1970a, b) at 100 to 400 MeV, and by Raisbeck and 

Yiou (19 75) above- 1 GeV. 

' 26 

The less abundant (11%) isotope Mg is a- significant source 
of gamma ray line emission at 1.809 Mev from deexcitation of 
the first level and by (p,, n) production of long-lived ^®Al 

Q ^ O 

which decays back to °Mg -through the T.809 MeV level- °A1 is 
also listed in Table 3- as an important long-lived gamma fay produc- 
ing radioisotope. The cross section for direct excitation of the 
1.809 MeV level, measured at 17.5 MeV by Crawley and Garvey (1967) 

is about 35% of that for excitation of the 1.369 MeV level in 

2 4 26 

Mg. In our siabsequent calculations we assume -that -the Mg 

(P/P' gQg) and^^Mg(a,a’'Y.^ gpg) ^^Mg cross sections are 

one half of ' the ^^Mg(p,p'Y 3 ^ 2 gg)^^Mg-and ^'^Mg(a,a' 368^^^^^ 
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cross sections, respectively, shown in Figure 8. The cross 

section for the '(•p,n) reaction leading to long-'lived Al has 

been measured' from 8 to 14 MeV by Wong et al. (196 7) and from 

‘6 to 52 MeV by Furukawa et al, (1971) , and is shown in- Figure 9. 

We also include for j completeness in °A1 production the' 'a 1 (p,pn) ^” a1 

cross section measured, by Furukawa. et al. (1971) . The production 
26 28 

of Al from Si spallation is discussed below. 

2 8 

'The stronest’ deexcitation line in Si is at 1.779 MeV. The 

2 8 

solid curve in Figure 10 is the cross section for the reaction Si 
20 ."*1 779 

(p,p') Si * . This cross section was measured by Conzett (1957) 

at 12 MeV, Crawley and Garvey (1967) at 17.5 MeV, deSwiniarski 

et al. (1973) at 25.25 MeV and Horowitz et al. (1969) at 100 MeV. 

Foley et al. (196 2b) ‘-measured the cross-section for 1.779 MeV 

gamma rays at proton energies of 141 MeV. Their result is shown 

by the open square in Figure 10. As with ^^Mg, the reason for the 

high value of this measurement is the contribution of the cascades 

2 8 

from the numerous {> 40) bound states’ of Si which deexcited 

mainly through the 1. 779 MeV level. The dashed curve in 

Figure 10 is our estimate of the cross section for. the reaction 
28 2 8 * 

Si(p,p'"Y^ ^^g) Si based on preliminary gamma ray data of Dyer 
and Bqdansky (private . communication , 1977) up to proton energies of 
25 MeV and the measurements of Foley et al. (1962b) at higher 
energies. 

2 8 

Two Other strong lines from Si are those at 6.878 and 5.099^ 

28 *6 879 

MeV from deexcitation of Si * . The cross section for proton 

excitation of this level, measured from 12 to 15 MeV by Shotter, 
Fisher-and Scott (19 70,), at 17:5 MeV, Crawley and Garvey (1967) and 
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at 100 MeV by Horowitz' et al. (1969) , is roughly 17% of that for 

op 2 R 

excitation of Si{p,p'Y^ Si. Using the branching ratios 

of de Meijer et al. (1975) and assuming a 20% cascade contribution 
from higher levels we estimate the 6.878 and 5.099 MeV gamma-ray 
production cross sections to be, respectively, 13% and 6% of 
that for 1.779 MeV gamma-rays in Figure 10 above an effective 


proton threshold energy of 8 MeV. 

28 28*1 779 

The cross section for the reaction Si (a, a') Si 


shown by the dashed-dotted curve in Figure 10. This cross section 
has been measured by> Blatchley and Bent (1965) at 5.6 MeV/nuclebn, 


by Kokame et al. (1965) at 7.1 MeV nucleon and by Rebel ^ al . 
(1972) at 26 MeV/nucleon. At 4.5 MeV/nucleon we use scaled data 
of Corelli et al. (1959) for the reaction ^^S(a,a’)^^S 
At lower energies, we again extrapolate the data according to the 


expected behavior near the Coulomb barrier. The dashed- crossed 


curve in Figure 10 is our estimate for the cross section of 
2 8 2 8 

the reaction Si (a, a' ^^g) Si. In analogy with the (p,p') 

reaction, we take a factor of 2 for the ratio of the cross 

^sections of the reactions ^^Si(a,a*Y]^ 779^^^^^ ^^Si(a,a') 

28 *1 779 

Si"" * at energies well above the Coulomb barrier energy. 

28 *6 879 

The cross section for a particle excitation of Si 


which leads to 6.878 and 5.099 Mev gamma- rays has been measured 


at an a energy of 7.1 Mev/nucleon by Kokame et al. (1966). As 

2 8 

with (p,p') excitation, these measurements show that the ‘ Si 
(a,a')^^Si cross section is roughly half that for ^^Si(p,p') 

oQ *1 779 

Si ■ . Using the branching ratios of de Meijer et al.(1975), 

we estimate the 6.878 and 5.099' MeV gamma- ray production cross 
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section from (a, a'') to be, respectively, 13% and 6% of that for 
/ 

1.779 MeV gamma- ray production. 

Foley et al. (1962b) observed a rich spectrum of proton 

''2 8 

induced spallation lines from Si. The strongest such line 

24 ... 

is at 1.369 MeV from Mg. The cross section for thxs reaction 

is larger by about a factor of 2 than, the cross section for the 

' 16 

production of the 4.438 MeV line from 0 at this energy. Based 

on this measurement, we ass^ame that the cross section for the 
28 24 * 

reaction Si(p,xy^ 369^ is, at all energies, 2 times the 

X6 12 

cross section for the reaction 0(p,xy^ 438^ ^ shown in 

2 8 

Figure 3. Other spallation lines of Si are at 1.634 and 

20 27 

2.613 MeV from deexcitations in Ne, at 0.780' and 0.957 from Si, 

2 7 

and at 0.844 and 1.014 from Al. From the data of Foley et a-1. 
(1962b) the cross sections for producing these lines, at an in- 
cident proton energy of ~ 150 MeV, are, respectively, 65, 40, 

25, 20, 30 and 40% of the 1.369 MeV line cross section. We 
use these ratios in our subsequent calculations, and we assume 
that they remain the same at all incident proton energies. Since 

Zobel et al. (1968) have observed the 1.014 MeV line from direct 
27 

excitation of Al, we also include in our calculations the 
27 27 

reaction ‘Al(p,p'y^ 014^ which we take an average 

28 28 

cross section of 50% of the Si(p,p'y^ Si cross section, 

based on the Zobel et al.(1968) measurements. 

As with Ne and Mg, we take into account the many other weaker 
lines from Si and its spallation products by considering the 
total gamma-ray production cross section. There are no direct 
measurements of this cross section for Si. The diamonds in 
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Figure 10 are estimates obtained by interpolating the data of 
Zobel et al. (1968) , -as was done for He and Mg. The crossed' 
and crossed- dotted curves in Figure 10 show the total gamma-ray 
production- cross section in Si for proton and a particle induced 
reactions , respectively . 

22 26 

An additional contribution to Na and Al production is 

obtained from the spallation of Si. Measurements of the Si 

spallation production cross section for these two isotopes by 

protons^ of energy up to 52 MeV are reported and reviewed by 

22 

Furukawa et al. C1971). Higher energy measurements of Na 
production have been made at 50 to 150 MeV by Bintoot and Gauvin (1971) , 
at 100 to ,400 Mev by Korteling and Caretto (1970a, b), and above 1 
GeV by Raisbeck and Yiou (1975). These cross sections are shown 
in Figure 9. 

32 

The principal gamma ray line from deexcitation of S is at 

2.230 MeV from the first level. The cross sections for excitation 

of this level has been measured at 5 to 5.5 MeV by Oda et al. 

(1959) , at 5 to 14 MeV by Berinde, Neamu and Vladuca (1971) at 

7.6 to 14 MeV by Oda et al. (1960), at 17.5 MeV by Crawley and 

Garvey (1967), and at 155 MeV by Willis et al. (1968). These 

cross sections are roughly 70% of that for excitation of the 

2 8 

first level in Si at 1.779 MeV. Assuming that cascade con- 
tributions are comparable, we use cross sections for both 

^^S(p,p’Y 2 230^^^^ ^^S(a,a'.'Y 2 230^^^® ^ times 

28'' 28 28 28 
those for Si(p,p'Y 2 779 ^ Si(a,a*YjL 779 ) Si, respectively. 

32 

Spallation of S also leads to strong gamma ray lines at 

9 fi * 

1.779 Mev from S(p,pa) Si and at 1.249, 1.266, 2.029, 2.034, 
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•30 * ^o -:j-| * 

2.232 and 2.234 MeV from '^S(p,2p) P and (p,pn) . Foley 
et al. (1962b) have measured -xcross sections for the 1. 779 Mev 
line and. for each of the pairs at ~1.2, ~2.0 and ~2.2 MeV and 

find them to be 160%, 330%, 120% and 140% of that for the 4.438 

< 

16 - 12 * 

MeV line from 0(p,px) C , respectively. For the latter pair 
we have subtracted the contribution of (pp') at 2.230 MeV measured 
also by Foley et al. (1962b) . Assuming similar excitation -func- 
tions and equal yields for each of the paired lines, we take 

•the cross sections for the 1.249 , 1.266, 1. 779, 2,029 , 2.034, 

2.232 and 2.234 MeV lines from spallation to be 1.7, 1.7, 

1.6, 0.6, 0.6, 0.7 and 0.7 times that for the 4.44 MeV line from 

16 ' 

O spallation (Figure 3) , respectively. 

32- 

For the unresolved line contribution from S we assume that 

for both proton and a particle excitations the cross sections 

2 8 

are the same as for Si given by -the crossed and crossed-dot-ted 
curves in Figure 7. 

40 

The strongest deexcitation line from Ca is at 3.736 MeV from 

\ 

the first gamma decay level. The cross section for proton excita- 
tion of this level has been measured at 14.6 and 17.3 MeV by 
Gray, Kenefick and Kraushaar (1965) , at 30 MeV by Ridley and 
Turner (see Satchler, 1973) , at 55 MeV by Yagi et al. (1964) , and 

at -155 MeV by Willis et al. (1968) and Roos and Wall (1965). 

< 

These cross sections lie essentially halfway between those for 

proton excitation of the first levels in ^^Si and ^^Fe or 150% of 
2 8 

•that for Si. Assuming that the cascade contribution is also 

2 8 56 

an average between -that in Si and Fe, we take the cross sections 
for '^^Ca(p,p'Y 2 ^*^Ca(a,a'Y 2 equal to 

^^Si (p,p' Yj^ 779 ) ^^Si and the ^^Si(a,a'Yi 779 ^^^®^ 


1.5 times the 
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cross sections, respectively. For the unresolved lines of Ca 

we again assume a cross section that is 1.5 times the corresponding 

. 28 ^. 

one in Si. 

•56 . 

■The strongest deexcitation line of Fe is at 0.847 MeV. The 
cross section for direct excitation of -bhis level by the reaction 
^^Fe(p,p') ^^Fe 0*847 shown by the solid curve in' Figure 11. 
From 3.5 to 6 MeV the data is from Nichols et al. (1969), at 6 MeV 
from Andronov et al. (1970, at 11 MeV from Perey et al. (1970) 
at 19.6 MeV from Hendrie et al. (1969), at 17.5 MeV from Peterson 
(1969), at 30.3 MeV from Karban et al. (1970, and at 49.35 MeV 
by Mani {19 71) . 

The cross section for -bhe reaction ^^Fe(p,p’Y_ o/n)^^^’© was 
measured by Zobel et al. (1968) at 15.7 and 31.4 MeV, and by 
Chang et al. (1974) and Jastrzebski et al. (1976) at 100 MeV. 

Their data are shown in Figure 11 by open squares and triangle 
and circle. As can be seen, there is an inconsistency of a 
factor of 3 between the data of Chang et al. (19 74) and Jastrzebski 
et al. (1976), which hopefully will be resolved by new measure- 
ments (N. S. Wall, private communication 1977) . 

The cross sections for producing 0.847 MeV gamma rays shown 
in Figure 11 are much larger than the cross section for the reaction 
Fe(p,p') Pe * due to -the contribution of cascades from 
higher excited levels of ^^Fe. This nucleus has more than 100 
bound states (Nuclear Data Group 1973) which are known to decay 
almost exclusively by cascades via the 0.847-MeV level. From 
Peterson (1969) at 17.5 MeV, and Mani (1971) at 49.35 MeV, we 
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see that the cross section for the excitation of each excited 

state is on the average about 10% of the. cross section for the 

excitation of the 0.847 MeV level. Thus at these energies the 
<* 

cross section for 0.847 MeV gamma ray production should be 

56 56 *0 847 

larger by about an order of magnitude than the Fe(p,p*) Fe 
cross section. In our calculations we use the cross sections 
given by the dashed curve in Figure 11. The comparison of this 
curve with the preliminary data up to 32 MeV based on gamma- ray 
measurements (P. Dyer, and D. Bodansky, private communication 1977) 
gives a good agreement. 

The cross section for producing 0.847 MeV gamma rays from ^^Fe 
5 6 

"^0.847^ Fe has been measured by Zobel et al. (196 8) who 
find a value of 79Q+ 430 mb at 14.3 MeV/nucleon. Because of the effei 
of the Coulomb barrier at low energies', we assume that this cross 
section is proportional to the dashed curve in Figure 11, with 
a constant of proportionality of 1,5 based" on"the "data "o^ Zobei. 
et al. (1968) . 

56 

Another strong line from Fe is at 1.238 MeV produced by 

transitions from the 2.085 MeV to the 0.847 MeV levels; these 

transitions are always followed by the emission of 0.847 MeV 

; 

photons. In (p,p') reactions the strength of the 1.238 MeV line 
relative to that of the 0.847-MeV line is 57% at 15.7 MeV (Zobel 
et al. 1968) and 45%. at 100 MeV (Chang et al. 1974). At 31,4 
MeV Zobel et al. (196 8) give a cross section for the 1'. 238-MeV 
line larger by about a factor of 2 than that for the 0.847 MeV 
line. This is inconsistent with the fact that every 1.238 MeV 
photon should be followed by a 0.847 MeV photon, but because' of 
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poor detector resolution, the measureraerits of Zobel et al. (1968) 
of ~l.'3-MeV photons probably included a large contribution from 
spallation lines at energies other than 1.238 MeV. Since spalla- 
tion reactions contribute less at lower energies, the measure- 
ments at 15.7 MeV probably give a more correct value for the 
1.238-MeV line cross section. At 100 MeV, the measurements of 
Chang et al. (i974) have adequate energy resolution to resolve 
the 1.238-MeV line from spallation lines which are quite dominant 
at this ener^. In the subsequent calculations we have taken a 
value of 0.5 for the ratio of 1.238 to 0.847 MeV gamma-ray pro- 
duction. This value is roughly consistent with the preliminary 
data of Dyer and Bodanksy (private communication 1977) . 

In addition to emission at 0.847 and 1.238 MeV, deexcitations 
of excited states in lead to gamma ray production also at 

a variety of other energies (Table 2) . We have included in 
Tables 1 and 2 all lines whose intensities are^ expected to be 
more than 10% of the 0.847-MeV line intensity. From Rao (1970) 
and Chang (1974) such lines are at 1.772, 1. 811, 2.094 and 

2.113 MeV, with intensities relative to that of the 0.847-MeV lijie 
of approximately 14, 50, 12 and 10%, respectively. As with the 
1.238-MeV line, we assume that -the energy dependences of the 
cross sections of these lines are the same as that of the line at 
0.847 MeV. 

5 6 

Spallation of Fe leads to several important prompt gamma ray lines 
(Table 2). At low energies (<20 MeV) , Dyer and Bodansky (private 
communication ,1977) find a strong line at 0.812 MeV from the 
reaction ^®Fe (p,n) ^^Co whose" cross section at 10 MeV 
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is about 30% of tliat of the reaction ^^Fe (p,n) ^^Co given in 
Figure 12.' We assume -this ratio at all energies in our subsequent 
calculations. The line at 0.812 MeV implies an additional de- 
excitation leading to 0.15 8 Me'V photons, whose intensity is 
at least equal to that of the 0.812 MeV line. 

At 100 MeV Chang et al. (1974) find the strongest spallation 

56 52 * 

line is at -1.434 MeV from the reactions Fe(p,x) Cr and 
56 52 * + 52 

Fe(p,x) Mn (e ;e) Cr. The latter reaction populates the 
. . 52 

xsomeric state of Mn which decays directly into the excited 
52 

state of Cr. There' are no measurements of the cross section 
52 

for Cr* production at other energies. However, we can estimate 

16 12 

this cross section by comparing it with the reaction 0(p,xY^ C 

in Figure 3. By using the results of Jastrzebski et al. (1976) , 

5 6 5 2 

we take the cross section for the reaction °Fe(p,xy^ 434^^ 
equal to 3 times the values given by curve 1 in Figure 3. 

The other prompt spallation lines of ^^Fe shown in Tables 1 
and 2 are those whose cross sectionsat 100 MeV.- are at least 10%' 
of the 1.434 MeV line cross section. From Chang et al. (1974) 
such lines are at 0.092, 0.412, 0.477, 0.744, 0.931, 

1.312, 1.334, 1.370 and 1'. 40 8 MeV with cross sections of 15, 15, 

15, 20, -65, 30 , 15, 10 and 60% of the 1.434 MeV cross section. 

We use these ratios in our calculations at all energies. 

As with Ne, Mg, and Si, the total gamma-ray production cross 
section from ^^Fe is considerably larger than the cross section 
for the excitation of any single line from this nucleus. In 
Figure 11 the closed squares represent data for this total cross 
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section (2obel et al. 1968) . The diamond is based on an extrapol- 

X6 2 7 

ation of the measurements of Zobel et al. (1968) for O and Al 
at 56 MeV. We have assumed isotropic gamma-ray emission in trans- 
formi-ng the differential cross sections given by these authors 

I 

into total cross sections. The crossed line in Figure 11 is 

our estimate of the total gamma-ray production cross section by 

proton- induced reactions on °Fe. For the unresolved gamma- ray 

. 56 

emxssion from a- particle- induced reactions on Fe, we take 
cross sections that at all incident energies are 3 times larger 
than the crossed line in Figure 11, based on a measurement of 
Zobel et al. at 14 MeV/nucleon. 

Spallation reactions on ^^Pe also produce a rich spectriim of gaitma ray 

lines from the decay, of long-lived radioisotopes. The most 

. 56 55 54 52 

important such isotopes are Co, Co, Mn, and Mn, 

along with ^Be, 22jjg^ have already been discussed above. 

56 I 

Decay of Co leads to strong line emission at 0.847 and 1.238 

MeV with weaker lines at 2.599 , 1.772' and 1.038; ^^Co decay gives 

a strong line at 0.9 31 MeV and somewhat weaker lines' at 1.40 8 
54 

and 0.477' MeV; Mn decay yields just a single strong line 
52 

at 0.835; and Mn decay produces three strong lines at 1.434, 

0.936 and 0.744 MeV. These lines, included in Table 1, have 

branching ratios of at least 10% (Bowman and MacMurdo 1974). 

The cross sections for production of these isotopes from 

56 

proton spallation of Fe are .shown in Figure 12 . Measurements 
of these and other ^^e spallation cross sections have been sum- 
marized and discussed by Bradzinski et al. (1971) for proton 
energies greater than 10 MeV. Measurements of the cross sections for 
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56 

Co production at lower energies have been sToramarized by Jenkins 
and Wain (19 70) . 

Gamma ray line emission from decay of other relatively long- 

51 48 46 

lived isotopes, such' as Cr, V and Sc, is not significant 
compared to those above,' either because the peak production 
cross section is < 10 mb' or' the gamma ray emission probability dur- 
ing decay. is too small (< 10%). 
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III. GA^MA. RAY SPECTRA AMD INTENSITIES 

In this section we discuss the shapes of the gamma ray lines 

and we carry out calculations of their intensities for various assump- 
tions on the spectrum and conqjosition of the energetic particles. We 

also -treat the formation of gamma ray lines in interstellar grains. 

\ 

In Ilia we discuss in detail the shape of the 4 .438 MeV line 
from (p,p^) reaction on V. Detailed measurements of the' spectrum 
of this line have recently been made in the laboratory and we compare 
our calculations with these measurements. In Illbwe treat the effects 
of interstellar grains on gamma ray line widths, and inlUcwe discuss 
the kinematics of the various other lines considered in .this paper. 

In Illd we present numerical calculations of gamma ray spectra and in- 
tensities from energetic particle interactions for a variety of spectral 
parameters and compositions. 
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a. Profiles of the 4.438 and 6.129 MeV Lines from (p,p*) Reactions 

Consider the 4.438 MeV gamma rays resulting from the interaction of 
12 

a proton with a . C target. We take the z axis in the direction of the 
Incident proton. The probability of photon emission per second into 

solid angle dcosBgd^Q from interactions in which a proton of energy 

12 *4.439 * 

E produces a C nucleus with recoil velocity in dcos0^d<^^ in the 

center-of-mass frame of the reaction is given by 

= nj.v ^fcCE,e*)dcos0*d(t)j. -g(E,0*,0Q,(})^-(j)^)dcos9^^d<j)Q. (1) 

Here n^. is the number density of the carbon- target, v(E) is the velocity 
of the proton, da/dfi*'is the center-of-mass differential cross section, 
and g is the angular distribution of the gamma rays which can depend on 
all the variables shown in equation (1) . Both 0* and 0^ are measured 

s 

with respect to the z axis, while the azimuthal angles (j)^. and (f)^ are, as 
usual, given in the (x,y) plane. 

A 

The gamma ray energy, Ey, corresponding to the variables E, 0j-, 

®o ’ ^ ^ single valued function of these variables through the 

following set of equations 


Ey = E^/ (y^- , . 

^r = ’ 

cos0y = cos0^cos0^+sln0j.sin0QCOs((|)j--^^) , 

/ 1/2 

cosOj. = [y*(Yc"^) YcW* -1) cos0*]/(y^-1) , 

Vc ■ 

Y.* ' 

‘ r cm 0 4 o cm 

^‘cm = E(mi + m2)^ +2 m 2 E]^'^2^ 


( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


( 9 ) 
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12 

Here is the Lorentz factor of the center of mass of the p C 

* 12 
system, y^ and y_^ are the Lorentz factors of the recoil C in the 

cent er-of -mass and the laboratory frames, respectively, 0^ is the 

angle between the velocity vectors of the incident proton and the 

recoil in the laboratory frame, 0^ is the angle between the velocity 

vector of the recoil G and the direction of the gamma ray in the 

laboratory frame, m^^ and m^ are equal to the proton mass mp, m^ is the 
X2 

mass of the • C, and m^ = m 2 + Ae where As = 4.439 MeV. Because of 

12 

additional recoil of the C during photon emission, the photon 

12 *4.439 

energy in the rest frame of the G . is given by 

= AS[1- Ae/(2m3C^)]. 

12 / 12 *4.439 

The total cross section for the reaction C(p,p ) G -* was 

discussed in detail in Section IIb» For the differential cross section 

of this reaction, we use the data of Peele (1957) at E = 14 MeV 

for the proton energy range 12£E<20 MeV, that of Dickens, Haner and 

Waddell (1963) at 24.1 MeV for 20£E<35 MeV, the data of Fannon. et al. 

(1967) at 50 MeV for 35^E<70 MeV, and that of Strauch and Titus (1956) 

at 96 MeV for ,E^70 MeV. Below 12 MeV we assume that da/dG^- is isotropic 

For the Monte Garlo simulation that we use for the evaluation of 

gamma ray spectra we shall need numerical values for the probability 

that cos0^ is less than a given value. Using the above data on da/dO* 

we show these probabilities by the solid curves in Figure 13. ,As can be 

12 * 

seen with increasing proton energy the G nuclei tend to recoil more 
towards the backward direction in the center-of-mass frame. For example 
half the nuclei have recoil angles, 0*, greater than 90° if E<3.2 MeV, 
whereas if E = 96 MeV half have 0*>155°. 

We proceed now to discuss our assumptions for. the gamma-ray angular 



39 


distributions g. In previous treatments (Bamaty and Crannell 1976, 

R^naty, Kozlovsky and Suril977) it was assumed that the photon angular 

12 *4 439 

distribution is isotropic in the rest frame of the C ' nucleus. 

> . 

This assumption, however, cannot explain the profiles of the 4.438 MeV 

line obtained in accelerator' experiments in which a high resolution 

gamma ray detector was placed at a fixed angle to the direction of the beam. 

For example when this angle is 90®, the line profile has a. deep minimum 

at 4.438 MeV, and two symmetrically located maxima whose separation, AKy, 

depends on the energy of the protons: AE^ 90 keV for E = 50 MeV (N. S. 

Wall, private connaunication 1977) and AEy 60 keV for both E = 23 MeV 

(Kolata, Auble and Galonsky 1967) and E = 16 MeV (P . Dyer and D. Bodansky, 

private communication 1977). Below 10 MeV the two peaks can no longer 

be separated (P. Dyer and D. Bodansky, private communication 1977). The 

ratio between the maximum and minimum intensities is about a factor of 

2 for all proton energies greater than 16 MeV, 

Following Kolata et al. (1967) , we can account for this splitting 

effect by assuming that the magnetic sublevels of the 2"^ state of the 

12 12 12 439 

C nucleus (at 4.439 MeV) are unequally populated in the C(p,'p') C 

reaction. In this case the angular distribution of the gamma-ray emission 

can be written as 

g = I * 

Here, cos0^ = slnG^ sinCcj)^.-^^) , where 0^ is the angle between the 

gamma ray and the normal to the reaction plane defined by the velocity 

12 * ' 

of the proton and the recoil C nucleus , and the are proportional 

I 

to vector spherical harmonics (e.g. Jackson 1962, page 551). The 

are the probabilities for populating the sublevels , and they can depend 
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on the proton energy E and the scattering angle 0 ^,. Schmidt et al. 
(1964) have measured these probabilities at, 10. 3 MeV. We have approxi- 
mated their results by the expressions 


Sq = 1 - (Si+S_3^) - (S2+S_2) (11) 

S^+S_j^ = 82 + 8^2 ~ 0*5exp[ (cos 0 *-l) /1.4] , (12) 

and we have assumed that these coefficients are independent of energy. 

Eqiiations (11) and (12) coupled with the angular dependence of 

do'/dO (Figure 13) lead to -the reduced population of the m = ^1 and +2 

sublevels, since and S ^2 largest at cos 0 ^ = 1 where da/d^ is the 

smallest. This effect can explain the dip at the center of the '4.438 MeV 

line, since photons at this energy result from reactions whose reaction 

plane is perpendicular to the direction of observations, and because 

✓ 

only the m = +1 terms in equation ( 10 ) produce radiation in the 0 ^ = 0 “ 
direction. 

We use a Monter-Carlo simulation to integrate equation (1) over 
cos^ and 4>j~ for- constant E, 0 q and ({)q. We choose a large nttmber of 
pairs of independent random numbers and R 2 , each of which are uniformly 

distributed from 0 to 1. For every such pair we solve for cos0j. and 
(pj. from the equations 


= [l/cr(E)] / ‘^‘"®®r(dCT/dff)dfi* 


(13)' 


^2 " ^r/2TT, 


(14)> 


where the right hand side of equation (13) is plotted in Figure 13. 

We then evaluate Ey from equations (2) through (9) and g from equations 
(10) through (12) . '■ The probability of observing gamma rays, of energies 
between Ey to Ey + AEy is then proportional to the sum of all the g’s 
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for -which is in the range. 

The results are shown hy the solid curves in Figure 14 for 0^ = 90° 
and three values of E. As can be seen, line splitting with strong dips 
does indeed follow from the numerical calculation; both the depths 
of the dips and the separation between the peaks are essentially 
consistent with the accelerator data of N. S. Wall (private communica- 
tion 1977) P. Dyer and D. Bodansfcy (private communication 1977) and 
Kolata et al. (1967), ' 

Also shown in Figure 14 (dashed curhes) is the line profile 

resulting from isotropic gamma-ray emission (S^ = S 2 ^+S_-jl = S 2 +S„ 2 ) ^.t 

E = 50 MeV and 0 q = 90°. As can be seen, this profile also shows a 

small amount of line splitting caused by the preferential emission of 
12 ^ 

C .nuclei into the backward direction in the center-of-mass frame. 

This splitting, however, is insufficient to account for the large 
splitting and deep dips observed in the accelerator data. 

If we now consider gamma ray emission from nuclear interactions of 
an isotropic distribution of energetic particles, then in addition to 
the integrations over 0 ^. and cfij-, equation ( 1 ) should be integrated over 
9 ^ and (j)Q. Furthermore, for a distribution of particle energies, an 
integration over E is also required. We carry out these integrations by 
choosing two additional random numbers, R 3 and R^, and by solving for 
E and cos 0 q from the. equations 

R 3 = vN (E'’)a(E')dE*' - (15) 

o ^ 

R 4 = (l+cos 0 o)/ 2 , (16) 


where N^(E) is the number of energetic protons per unit energy, and 
C is a normalization constant; there is no need to evaluate <|>o explicitly 
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because equation (1) depends only on We then determine the gamma- 

ray spectrum as before when E and 0^ were kept constant. 

The results are shown in Figure 15 for an energy spectrim Np(E) 

2 ' " ^ 

E and two choices of the angular distribution of the gamma rays . 

The solid line is for isotropic gamma ray emission, and the dashed line 
is for the Sm's as given in equation (11) and (12) . As can be seen, 
there is no substantial difference between these two spectra. In 
particular the full widths at half maximum (FWHM) of these spectra are 
essentially identical. This result implies that the effects of the 
anisotropic gamma-ray emission which are very noticeable in the case of 
a proton beam, are essentially smeared out when the distribution of the 
energetic particles is isotropic. Therefore, in the subsequent discussion 
and calculations we consider only isotropic incident particle distributions 
and assume isotropic gamma— ray emission for the 4.438 MeV line as well 
as for all the other lines discussed in Section II- Some other observable 
features in gamma-ray line shapes resulting from anisotropic .particle — 
distribution in solar flares have been discussed previously (Ramaty 
and Crannell 1976, Kozlovsky and Ramaty 1977) . 

For the calculation of the 6.129 MeV line profile in (p,p^) reactions 
we use data, for dcr/dQ* measured by Crawley and Garvey (1967) at 17.5 MeV, 
and by Austin et al. (1971) at 17,’ 29.5 and 46.1 MeV. The integrals required 
in equation (13), based on these data, are shown by the dashed curves in 
Figure 13. As for the 4.438 MeV line, we assume that da/dfi is isotropic 
for E<L2 MeV, and we use the 17 MeV data for the range 12^E<20 MeV, the , 

t '* 

24.5 MeV data for 20^E<35 MeV, and the 46.1 data for 35^E<70 MeV. 

I 

12 

Because of the lack of data at higher energies, we have used the C data 
at 96 MeV (Stranch and Titus (1956) also for ^^0 in the range Es70 MeV. 
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b. Lines from Interstellar Graias 

In calculating the profiles) of the various gamma ray. lines discussed 
in II, we have to consider the effects of interstellar grains which can 
lead to very narrow line emission (Lingenfelter and Ramaty 1977)-. This 
line component is produced by deexcitation of toterstellar grain nuclei. 

As we have seen in 11a, the energy, E^, of a gamma ray emitted 
from a given nuclear level. is Doppler shifted from the rest energy 
due to the recoil energy, .Ej., acquired by the excited nucleus in the 
nuclear interaction. However, if this interaction takes place in a grain, 
at the time of photon emission the recoil energy, E^, .could be smaller than 

I 

its initial value, E^, because of the energy lost by the excited nucleus 

in the solid prior to deexcitation. The width of the gamma ray line could, 

therefore, be substantially decreased provided that the mean life of the 

level or of its radioactive parent (Table 1) , and the distance from the 

site of the interaction to the grain edge along the direction of the 

recoil, .are longer, respectively, than the slowing down time and the 

stopping range of the excited nucleus in the grain material. 

The rate of energy loss, dE^/dx, the mean rectilinear stopping range, 

•<X(Ej.)>, and the fluctuation in X due to straggling, (AX’/X) are shown in 

Figure 16 for ^^0, and nuclei slowing down in water assuming 

that the grains are predominantly ice. These quantities, however, are 

not strongly dependent on the assumed grain composition. These curves 

were derived by Bussard (1978) from the compilations of Northcliffe and 

Schilling (1970) and Winterbon (1975) . Typical recoil energies from 

5 6 16 

(p,p^) range from about a keV for °Fe to several tens of keV for 0. 

-3 

For a density of 1.5g cm such nuclei slow down in less than about 

-13 —4 

7x10 sec over distances less than <~10 cm. By comparing with the mean 
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lives listed in Table 1 and the expected sizes of interstellar grains, it 
is evident (Lingenfelter and Ramaty 1977) that many lines could have 
significant very narrow components. 

In our calculations we consider spherical grains having an exponential 

% 

size distribution and containing half of the interstellar C, N, 0, Ne and 
S and all of the Al, Mg, Si, Ca, and Fe. For such gradLns, the probability 
that a nuclear interaction takes place in a grain of radius between a and 
a+da is • 

P(a)da = (6a^^) a^exp(-a/aQ)da. (17) 

In our calculations we take a^ = SKlO^^cm. We also asstime that the 
, energetic particle flux is isotropic, and that the matter density and 
composition are constant throughout the grain. Then the probability that the 
recoil nucleus traverses a distance to the edge of the grain between r and 
r+dr is 

P(r)dr = (12a2-3r^)/(16a^)dr. ' (18) 

In the Monte Carlo simulation for the calculation of the total 
gamma-ray spectrum we take into account the effects of interstellar grains 

for all lines whose nuclear levels or radioactive parents have mean lives 

-13 ' 

longer than 10 sec (Table 1) . Each event in the simulation corresponds 

to a nuclear interaction in some target nucleus which has a recoil energy 

Ej., calculated from Equation (‘4,). If this interaction takes place in the 

grain, we evaluate a grain radius from Equation (17) , a distance to the 

grain edge from Equation (18) , a stopping range X = <X> + Ax from Fi^re 16, 

"I 

a AX from a gaussian distribution, <and a lifetime, t^, from an exponential 
distribution with mean life given in Table 1. If the nucleus stops' in the 
grain (r>x) , E^ is evaluated from 
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r^''(pv)"W(dEj./dx) = t., (19) 

e;. 

where p is the density and v the recoil velocity; this case includes 
the possibility that = 0, i.e. the nucleus stops before it deexcites . 
Equation (19) is. also used when r < X, provided that the nucleus deexcites 
before it reaches the the grain edge, i.e. when 

(pv)~^dE/(dE^/;dx) > t^. (20) 

E^^edge 

Here, E^, ® , the energy of the nucleus at the edge, is obtained by solving 

^<X(Ej.®^®®)> = <S(E^)> -r. (21) 

I If condition (20) is not satisfied, then E^' = i.e. the nucleus 

leaves the grain in an excited state and deexcites later in the inter- 
stellar medium where its energy remains essentially An exception 

26 3.3 

is A.L, whose mean life, 3.4x10 sec, is longer than the stopping time 

of this nucleus in the interstellar medium. Thus, for the 1.809 MeV line 

26 ■' 
from ki, decay we take = 0, independent of the grain size. 

Having evaluated the effect of the grains on the recoil energy of 

the excited nucleus, we evaluate the gamma ray energy, E^, from 

Equations (2) through (9), with the transformation E^-»E^' as described- 

above. 

c. ’ Kinematics of Other Nuclear Lines 

We proceed now to evaluate the spectrum 'that results from the 
superposition, of all the gamma ray emission produced by the reactions 
discussed in Section II. 

We first discuss reactions induced by energetic protons and alpha 
particles. The data and assumptions used to calculate the profiles of the 
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4.438 and 6 ,129 MeV line from (p,p’) reactions were given in (Ilia). 

For all other (p,p^) and (p,n) reactions we assume that the angular 

12 ^ / 12 * 

dependence of da/dQ*' is the same as that for .the C(-p,p •) C reaction 

shown in Figure 13. For (a,or^) reactions, we assume that du/d^ is 

isotropic. An can he seen from the cross sections of II, (e.g. Figure 2 

for the reaction ^^C(q;,c/ )^^C*) , (.a,a' ) reactions are important only at 

projectiles energies below about 10 MeV/nucleon, In this energy range, 

the data of Mitchel, Carter and Davies (1970) , Morgan and Hobie (1970) , 

Blatchley and Bent (1965) and Correlli, Bleuler and Tendam (1959) 

indicate that the (q;,q'^ ) differential cross sections have considerable 

structure with several maxima and minima, but on the average they can be 

approximated by an isotropic distribution. As discussed by Kozlovsky and 

Ramaty (1977), it can be assumed that da/dif is isotropic for the reactions 

^He(or,n)^Be* and ^He(af,p)^Li*. 

To calculate line profiles from spallation reactions, we in 

principle need data both on da/dQ Eand on the distribution of Ys » 

because in reactions with more than two bodies in the final state, this 

quantity is not a unique function of the incident energy, E. But in the 

absence of data on this distribution, we have assumed for multibody 

processes that y ~ Yc» i-e. that the velocity of the recoil nucleus is 

the same as that of the center of mass. This assumption is different 

from that made by Ramaty et al. (1977), and it leads to a larger width 

than found by these authors. From gamma-ray measurements, P. Dyer and 

D. Bodansky (private communication 1977) have found that the FWHM of the 

X6 

4.438 MeV line from 0 spallation is somewhat larger than 100 keV, a 
result which is more consistent with the present assunption. 

-13 

For all lines from levels with mean lives longer than 10 sec we 

44 
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carry .out the transformation E -»E ' to take into account the effects 

of interstellar grains as discussed in Illb. 

We estimate the energy distribution of the unresolved gamma-ray 

lines of Ne, Hg, Si and Fe from the measurements of Zobel et al'o (1968). 

These authors have measured the energy distribution of gamma rays from 

27 56 

the bombardment of A-L and Fe by protons of various energies. Their 
results are quoted by Shima and Alsmiller (1970) . These distributions 
peak between about 1 to 2 MeV, with the magnitude of maximum becoming 
more pronounced with increasing target atomic number and proton energy. 

As we have discussed in lie, the gamma rays at these energies arise 
mainly from transitions in high-lying levels of the nuclei and their 
spallation products. Indeed, from the measurements of Degnan et al. 

(1973) we see that for a variety of target nuclei the photon multiplicity 
Increases with increasing excitation energy of the residual nuclei. From 
these multiplicities it follows that, on the average, the photon energies 
are in the 1 to 2 MeV range. 

The energy distribution of the unresolved gamma-ray lines are 

shown in Figure 17 normalized to unit integral. The closed circles are 
56 

data for Fe based on the measurements of Zobel et al. (19 68)' at proton 

energies of 31.4 MeV, and the open squares show~the assumed gamma-ray 
28 27 

distribution for Si based on measurements of A^, at 30 MeV. We 
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A-t, data at lower energies for, Ne and Mg is based on the trend in the 
data of Zobel et al. (1968) according to which the peak in the 
1 to 2 MeV region becomes less pronounced for lighter targets and lower 
incident energies 

For the reactions induced by energetic heavy nuclei we proceed 
as follows. For the two-body reactions we use, with the appropriate 
transformations, the same kinematics and data as for the corresponding 
reactions induced by protons or alpha particles. For spallation reactions, 
we assume that the excited nuclexis moves in the same direction and 
with the same velocity as the incident heavy nucleus. For the unresolved, 
lines we use the same data as for the corresponding reactions induced by 
protons and alpha particles, but broaden the spectra shown in Figure 17 
by assuming, as we just did for the spallation reactions, that, the excited 
nuclei move in the same direction and with the same velocity as the 
incident nuclei. 

t 

-Interstellar grains have no effect on gamma ray lines from excited 
energetic heavy nuclei since the stopping ranges of ithese particles are 

much longer than the dimensions of the grains. The 1.809 "MeV line from 

26 26 27 28 

A-t decay produced by energetic Mg, A.t and Si nuclei, however, will 

26 

be very narrow, because the energetic (<100 MeV/nucleon) A-t stops" in 

-3 

the interstellar medium of density ~ 1 Hem in a time shorter than its 
mean lifetime. 

d. Nuclear Gamma Ray Spectra and Intensities 

Using the above data and assumptions, we have evaluated the total 
gamma-ray spectrum by en 5 )loying for each of the reactions given in 
Section II a Monte Carlo simulation sim il ar to the one described in 
Ilia for the '4.438 MeV line. 
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For the ambient medium we use the abundances of Ross and Aller (1976) 
shown in Table- 2. For the energetic particles, we assume that all 
particle species have the same energy spectrum, 

Ni(E) 5 , (22) 

Const ; E^Ec 

where N^(E) is the number density of species i per unit energy per 
nucleon interval around E, and s and E^, are spectral parameters; steep 
spectra correspond to large s and low E^, and flat spectra to low s or 
large E^. We use two sets of abundance ratios fdr the energetic particles 
at the same E: solar abundances from Ross and Aller (1976) (Table 3, 

column 1) and cosmic' ray source abundances from Shapiro, Silberberg 
and Tsao (1975)' (Table 3, column 2) . _ 

The results are shown in Figures 18, 19 and 20 for E^ = 0 and s = 2, 

3 and 4 and solar abtindances, and in Figures 21 and 22 for s = 3 and' 

E^ = 0 and 20 MeV/nucleon and cosmic ray abundances. There are 10^ 
photons in each of these spectra, and they are binned into energy 
intervals of widths ranging from 2 to -5 keV, as indicated in the figures. 
Such widths are consistent with the energy resolution of Ge detectors 
used in current and planned gamma ray astrophysical spectroscopy missions 
(Lingenfelter and Ramaty 1978) . 

For comparison, in Figures 23 and 24 we show gamma ray spectra, 

obtained by binning the photons into 100 keV intervals. In these figures, 

we also use cosmic ray abundances for the energetic particles, and s = 3 
and E^ = 0 or 20 MeV/nucleon. The solid lines show tota'l spectra, while 
the dashed lines are spectra produced by only the energetic heavy nuclei 
(Aal2) . 

Three line components contribute to the spectra of Figures 18 through 
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24. There is a broad con^onent from the deexcitation of energetic ' 
heavy nuclei which interact with ambient H and He, a narrow component 
from ambient nuclei excited by energetic protons and alpha particles, 
and a very narrow component from deexcitation of those excited .nuclei 
that can come essentially to rest before deexcitation because. the 
lifetime of the level, or its radioactive parent, is long enough and the 
ambient density high enough. The very narrow component is mostly due to 
the interactions of protons and alpha particles with heavy nuclei is 
interstellar grains (Lingenfelter and Ramaty 1977). 

The broad lines tend to overlap, and therefore only a few well 
defined features can be distinguished in this component. Such features 
are most noticeable for steep particle spectra and cosmic ray abundances, 
and can be seen, for example, in Figures 21 and 23 at about 0.85, 1.3, 

1.7, 4.4 and 6.2 MeV. 

Many narrow and very narrow lines are superiu 5 >osed on the continuum 
defined by the broad component. These are particularly numerous if the 
energetic particles have the same composition as the ambient medium 
(Figures 18, 19 and 20); if the relative abundances of heavy nuclei are 
enhanced as in cosmic ray sources, then the underlying continuum becomes 
higher, and therefore fewer narrow and very narrow lines can be seen 
(Figures 21 and 22) . 

The strongest narrow line is at 4.44 MeV resulting from deexcitations 
12 11 

of excited states in C and B (Table 1) . Because of the short mean 
lives of these levels, the 4.44 MeV line has no veTry narrow conqjonent. 
However, -the second strongest line, at 6.129 MeV, has a prominent very 
narrow component which can be clearly seen in Figures 18 through 22 . 

In addition, there are many other lines which exhibit very narrow components 
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The strongest among these, at 0.847 and 1.238 MeV from Fe, at 1.369 MeV 

i 

from ^^Mg, at 1.634 MeV from ^^Ne, and at 1.779 MeV from ^®Si, are 

visible even for a cosmic ray composition, as can be seen in Figures 

26 

21 and 22. The 1.809 MeV line from Af, decay is very prominent in Figure 22. 

This line is produced in the interstellar medium which slows down the 

26 28 27 26 

excited A-t -nuclei Resulting from energetic Si, A-t and ■ Mg before 

26 

they decay into the 1.809 MeV level of Mg. 

The production rates, q, of the 4.44, 6.129 and 0.847 MeV lines 

and of the total nuclear radiation, normalized to the energy density, 

W, in energetic particles and IH atom .in -the ambient medium are shown 

in Figure 25. These results are based on solar abundances for both the 

\ 

ambient medium and the energetic particles. For the 4.44, 6.129 and 
0.847 MeV lines, the q's of Figure 25 are those resulting from proton 
and alpha Interactions, i.e. the narrow component of the 4.44 MeV 
line and the sums of the. narrow and very narrow components of the other 
two lines. For the total, q is the sxim of the .production rates of photons 
of all energies in all three line components . 

As can be seen from Figure 25, the q/W*s are maximal for a given 
s when is between about 10 to 30 MeV/nucleon. This is due to the 
fact that the nuclear cross sections peak in this energy region. 

The production rates, q, can also be normalized to the energy 

deposition rate, W, of the energetic particles in the ambient medium. 

» '■ 

The ratios q/W are shown in Figures 26 and 27 for a neutral and ionized 

ambient medium, respectively. As in Figure 25, the q*s are the production 

rates of the narrow 4.44 MeV and narrow plus very narrow 6.129 MeV line, 

and both the ambient 'medium and energetic particles have solar abundances. 

The energy loss fates in a neutral medium of 90 % H and 10% He are from 

Northcllffe and Shilling (1970) and Barkas and Berger (1964) . The energy 

loss rates in an ionized medium are from Ginsburg and Syrovatskli (1964) 
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calculated for a temperature of lO^K and electron density 1cm . As car> 
be seen, in a neutral medium q/W is larger -b,y about a factor of 4 than 
in an ionized one, because for a neutral gas W is smaller by this factor. 
By comparing Figures 26 and 27 with Figure 25, we see .that the q/W’-s 
tend- to peak at higher than q/W. This is due to the fact that the 
energy loss rates of fast particles increase very- rapidly as their 
energies decrease. 

The ratio of the production rates of the narrow plus very narrow 

6,129 MeV line to the narrow 4.44 MeV line is shown in the upper panel 

of Figure 28 as a function of s and E for solar abundances . As can 
• c . 

be seen, this ratio is quite insensitive to the spectral parameters, its 
value ranging from about 0.3 to 1. 

The production of nuclear gamma ray lines by energetic particle 
interactions is accompanied by positron production from the decay of 
radioactive positron emitters. Such positron production has been treated 
in considerable detail by Ramaty et al. (1975) . We have updated the 
cross sections used by these authors, and we have added several new 
positron emitters such as ^^0 and ^^Co. We shall present a detailed 

discussion of positron production elsewhere. In the lower panel of 
Figure 28 we show the ratio of positron production to the narrow 4.44 MeV 
line production as a ftmction of s and for solar abundances. We note 
that q^./q/|^ is not a strong function of the spectral parameters. None- 
theless, for extreme flat spectra, this ratio is larger by about a factor 
of 6 than for very steep spectra, reflecting the fact that positron 
production does not drop off at high energies as does the cross section 
for 4.44 MeV line production (Figure 2). We note that q^ in Figure 28 
does not include positron production from -n+ decay. 
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IV. AN ASTROPHYSICAL SOURCE; NUCLEAR GAMMA RA.YS FROM THE INTERSTELLAR 

MEDIUM 

Detectable nuclear gamma ray line emission may be produced- by energetic 
particle interactions in various sources. Nuclear gamma -ray production in 
solar flares bas been treated in detail by Ramaty et al. (1975, 1977) and 
other sources which could be obseirved in gamma-ray lines have been reviewed 
recently by Lingenfelter and Ramaty (1978). In this section we shall apply 
the theory developed in the previous two sections to the nuclear gamma ray 
line emission resulting from low energy (<100 MeV /nucleon) cosmic ray 
interactions with the interstellar gas and dust. We calculate the expected 
line intensities 'and we compare our results with previous estimates, as 
well as with reported observations of nuclear lines from the galaxy. We 
shall also compare the calculated nuclear line intensities with other diffuse 
lines expected from annihilation of positrons from various sources and from 
the decay of long-lived (>10^yrs) radioisotopes synthesized in novae and 
supernovae. 

Gamma ray line emission from cosmic ray interactions in the interstellar 
medium should constitute a spatially diffuse galactic source with maximum 
intensities in the galactic plane in the general direction of the galactic 
center, reflecting the distributions of matter and energetic particles along 
the line of sight. The intensity in any particular direction depends 
primarily on the distributions of the density and composition of interstellar 
matter, and the energy density, spectrum and composition of low energy cosmic 
rays. The intensity of the very narrow line component also depends on ‘the 
size, composition and distribution of the dust grains. 

Although these distributions are at best only poorly known, the line 
intensities can be calculated from the measured (Fichtel et al. 1975; Paul 
et al. 1978) intensity of high energy galactic gamma rays and the comparison 
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of the relative emissivities of nuclear lines and these gamma rays. A 
large fraction of the high energy gamma rays are believed to be due to 
the decay of tt** mesons produced by the interaction of high energy (>100 MeV/ 
nucleon) cosmic rays with interstellar matter. The nuclear line intensities 
are further constrained by x— ray and low energy gamma ray observations. The 
X— ray observations can set limits on the column depth in the galactic plane 
of 0 and heavier elements through studies of the absorption of radiation 

I 

from -x-ray sources (Ryter, Cesarsky and Audouze 1975); x-ray observations 
can also limit low energy cosmic ray fluxes through investigations of the 
nonthermal 6.8 keV iron line produced by charge exchange of energetic iron 
ions with interstellar gas (Bussard, Ramaty and Omidvar 1978) . A further 
limit on the in dividual nuclear line intensities is given by upper limits 
on the total gamma ray emission from the galactic disc in the MeV region 
(Gilman et al. 1978). We shall discuss these constraints below together with 
much more model depdndent constraints related to the ionization of the inter- 
stellar medium and Li production. But first we set up the foimiallsm for 

calculating the intensities of the nuclear lines. 

> 

Let and q^o be the emissivities per E atom of gamma rays in a nuclear 
line i and from tt” decay, respectively, as functions of galactocentric .radius, 
u). The fluxes that result from these emissivities, integrated over galactic 
latitude, b, at a given longitude, i,, can be written as 

I.(^) = l/4w f^°“db ^dr np(«))q (a,) (23) 

^ -10“ o 1 

and 

I oCt) = q o(0)/4tt ^dr n^(m)q o(o))/q o(0), (24) 

tt ilQO i ^ 

where r is the- distance from the observing point along the line of sight, 
n is the number density of protons (in both atomic and molecular form) 
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at 0 ), and is the local n° decay emissivity (i.e. at ,(jo=o)g = 10 kpc) . 

22 i 

The relationship between r» and m is u) = (o)g 4r ~2u)gr cos^) ; the limits 
of integration on h are consistent with presently used gamma ray line detectors. 
Let W^) be the energy density in low energy cosmic rays as a function 
of galactocentric radius, and W(0) the local energy density in' these 
cosmic rays. If we assume that the low and high energy cosmic rays have 
the same spatial distribution in the galaxy, and if t-heir .energy spectra 
are spatially independent, then 

W(oj)/W(0) = q^o(u3)/(lo(0). (25) 

From equations (23) , (24) and (25) we obtain that 

I.(L) = W(0) I o(L)/q o(0)<q.A^>, (26) 

-l- TT TT J- 

where <q^/lt> is the emissivity in line i per H atom per unit -energy 
density averaged along the line of sight . 

We shall use equation (26) to calculate the fluxes of the various 

3 

gamma ray lines. For W(9) we take the nominal value of leV/cm , roughly 
equal to the measured energy density of the high energy cosmic rays. The 
energy density in low energy cosmic rays outside the solar system is 
not known because solar modulation excludes these particles from the inner 
solar system (e.g. Goldstein", Fisk and Ramaty 1970). 

The second term in equation" (26) , I -o/q o follows from high energy 

, Tf TT 

gamma ray observations and theory. The high energy gamma ray spectrum 
has been measured by Paul et ai. (1978) from several directions in the 
galaxy. We shall use their data from the -general direction of the galactic 
center, integrated over about +10“ of galactic latitude from 355“ to 15“ of 
galactic longitude. This data is shoxm in Figure 29. .The general longitude 
dependence of the galactic high energy gamma rays was determined by 
Flchtel et al. (1975). For galactic longitudes between about 330“ to 30“ , 



56 


the high energy gamm ray intensity is constant to within a factor 
of 2. At gamma ray energies* less than about 200 MeV the spectrum can 
be fitted by a n^-decay gamma ray spectrum plus an apparent electron 
hremsstrahlung' component which we shall discuss below. But above 300 MeV 
the measured spectrum is consisteiut with that expected (e.g. Stecker 1976) 
purely from TT^-decay gamma rays (see Figure 29). Normalization, of the 
measured flux at these energies to the calculated TT°~decay spectrum gives 
^ o(^,)/q 0,(0) = 5x10 H atom/ (cm rad) for longitudes Ji near the galactic 

TT TT 

center. 

The emissivities per unit energy density, <q^/\}>, depend on the 
compositions of interstellar matter and low energy cosmic rays, and the 

• i 

energy spectrum of these cosmic rays . We first consider the composition of 
interstellar matter. - 

There may be evidence fot large radial gradients in the jr^ative abundance 
of elements in our galaxy, based on observations of HII emission nebulae 
in both our own and other spiral galaxies (e.g. Searle 1971, Shields 
1974, Smith 1975, Peimhert 1975, and Peimbert, Torres— Peimbert and Rayo 
1978). Such .composition gradients are expected from galactic evolutionary 
models (Tinsley and Larson 1978). Within a few kiloparsecs/ of the Sun, 

Peimbert et al. (1978) find radial gradients for 0 and N of d log (p/H)/d(u = 
0.13 jH). 04 kpc ^ and d log (N/H)/dy> = 0.23 +0.06 kpc Measurements 
(Burton 1976) of the longitude-velocity distribution of CO emission 
suggest a maximum in the molecular density in a region at galactocentric 
radius of 4 to 6 kpc. Since it is believed (Stecker 1976) that a major 
fraction of gamma rays from Tr° decay are produced in this region, we ' 
calculate for O/H and N/H at cju ““ 5kpc. From the local abundance 

gradients, these ratios should be larger than -solar system abundances by 
factors of ,>^5 and .J.5, respectively. The solar system abundances and the 
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assumed ambient medium abundances at 5kpc are shown in the first and 

/ 

third columns of Table 3. In the latter, the abundances of all elements 

heavier than He, except N, are larger by a factor of 5 than solar system 

abundances; the N abundance is larger by a. factor of 15. 

The assumptions regarding the enrichment of heavy elements in the 

interstellar medium are limited by the column density . of these elements 

as determined by interstellar absorption of x' rays, and- by the lower limit 

on the hydrogen column density obtained from 21cm observations. The column 

19 2 

density of 0 is about 7x10 cm (Tucker et al. 1973, Ryter, Cesarsky and 

22 2 

Audouze 1975), while that of H should be, greater than /sJ.0 cm , the value 

/ 

obtained from 21cm emission. These two column densities would allow a 

maximal O/H ratio larger by about a factor of 10 than the solar O/H. The 

/ 

enrichment factor of 5 used in Table 3, however, is more consistent with 
the presence of comparable amounts of atomic and molecular hydrogen along 
the line of sight. 

We next consider the composition of the low energy cosmic rays. Since 
this cosmic ray component has not been directly observed, we base our 
assumptions about its composition on information obtained from the high 
energy component. The composition of high energy cosmic rays at their 
sources has been calculated by several authors (e.g. Shapiro et al. 1975). 

From these calculations we assume a low energy cosmic ray composition as shown 
in the second column of Table 3. It is believed that the cosmic rays observed 
near Earth reflect galactic conditions in a region extending to a distance 
of not more thani^J. kpc from the solar system (e.g. Ramaty, Reames and - 
Lingenf elter 1970) . In the fourth column of Table 3 we present a possible 
cosmic ray source composition at 5kpc, where we assume that the abundances 
of nuclei heavier than He are enhanced relative to local cosmic ray source 
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abundances by the same factors as the ambient medium, i.e., by a factor 
of 5‘, except for N whose abundance is enhanced by a factor of 15. 

In our subsequent calculations we shall consider four combinations 
of abundance ratios, which we shall refer to as cases (1) through (4), For 
case (1) both the ambient medium and low energy cosmic rays have solar 
abundances; for case (2) the ambient medium has solar abundances and the 
cosmic rays have local cosmic ray source abundances; for case (3) both 
the ambient medium and the low energy cosmic rays have the same composition 
as the enhanced ambient medium composition at 5kpc; and for case (4) the 
ambient medixim has the same composition as for (3) , but the cosmic rays 
have an enhanced cosmic ray source compesition. 

The dependence of q^/W and q^^/W on the spectral parameters of the 
energetic particles, s and E^, was treated in Illd for solar abundances, and 
the results are shown in Figures 25, 26 and 27. In Table 4 we present 
<q/tt> and <q/^' for the narrow and broad 4.44 MeV lines for cases (1) through 
(4), and s = 4 and = 20 MeV /nucleon which essentially maximize these 
emissivities'; W in this table is calculated for a neutral ambient medium. As 
can be seen, q/w and q/W increase, in general, with increasing relative 
abundance of the heavier nuclei. But this increase is not linear, because 
for large heavy element abundances in the - energetic particles, a large 
fraction of W and W are due to the heavy nuclei themselves. The contribution 

of nuclei heavier than He to the total W is about 1.5%, 14%, 8% and 47% for 

/ * 

cases (1), (2), (3) and (4), respectively, and their contribution to W is 
7%, 51%, '29% and 84%, again for cases (1) through (4), respectively. Thus, 
for case (4) , for example, half the energy density and almost all the energy 
loss of the energetic particles is due to nuclei heavier than He. 
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Also shown in Table 4 is the ratio of positron production to the 
total (narrow + broad) 4.44 MeV line production. As can be seen, this 
ratio is essentially independent of the assumed abundance combination, 
since both the positron emitters and the 4.44 MeV photons are produced pre- 
dominantly from C, N and 0. 

For the , coii 5 )arison of the calculations with the data of Haymes et al. 

(1975) and Gilman et al. (1978), we also need the production rates of the 
total nuclear radiation, q(tot), and the production rates of gamma rays 
in the energy bands from 0.6 to 5.2 MeV, q(0..6-5.2), and 4.09 to 4.79 MeV 
q (4. 09-4. 79) . For s=4 and Eq = 20/MeV nucleon and' cases (3) and (4), respe( 
tively, <q(tot)/l^ = 1.4x10 and 6.0x10 photons/ (sec Hatom eV cm ^); 
q(0,6-5 .2)/q(tot) = 0.72 and 0.80; and q(4.09-4.79)/q (tot) = 0.17 and 0.085. 

Having discussed the three terms that are' required to evaluate equation 
(26) , we now proceed to calculate' the fluxes of gamma ray lines from the 

3 

general direction of the galactic center. By substituting W(0) = 1 eV/cm 

'20 2 

and I o(i)/q o(^) = 5x10 Hatom/ cm rad into this equation, we evaluate 
TT/ n 

the fluxes of the narrow 4.44 MeV line, the narrow plus very narrow 6.129 MeV 
and 0.847 MeV lines, and the 0.511 MeV line for cases (1) through (4). The 
resvilts are given in Table 5 for s=4 and = 20 MeV /nucleon. For the 
0.511 MeV line, we take a photon-to-positron ratio of 0.6. This ratio is based 

I 

on the calculations'- of Bussard, Ramaty and Drachman (1979) who find that in 

the interstellar medium about 93% of the .positrons should annihilate from 

bound states of positronium with a 0.511 MeV photon-to-positron ratio of 0.5,'" 

and the remainder by direct annihilation with a ratio of 2. 

The calculated line intensities given in Table 5 can be compared with 

earlier calculations. For the narrow 4.44 MeV line, for exan^le, Meneguzzi 

A 2 

and Reeves (1975) find an intensity of 5x10 photons/cm sec sr MeV if the local 

3 

local low energy cosmic rays have an energy density of 0.4 eV/cm . 
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Since the width of the 4.44 MeV line is 110 keV and the thickness in 

latitude of the galactic ridge is about 4°, this intensity yields a 
“6 2 

fliix 'of ^xlO photons/,(cin .sec rad) which- is- smaller by about a factor of 

-5 2 

3 than our predictions of (1.2 to 1.3)xl0 photons/cm sec rad for local 
conditions (cases 1 and 2). This difference is due to the lower energy 
density used by these authors, and the lower q/W resulting from their 
choice of spectral parameters (E^ = 5 MeV /nucleon and s = 3.5 and 5.5). 

Meneguzzi and Reeves (1975) chose these values to limit the rates of ioniza— 

7 A 7 

tion of interstellar H and ^Li production (by .the reactions ^He(o!jP) Li and 

4 7 -15 -1 -27 -1 

He^,p) Be) to 2x10 (H atom sec) and 3x10 (H atom sec) , respectively. 

For the choice of parameters of Table 4 (s=4, E^ = 20 MeV/nucleon) and W(0) = 

1 eV/cm^, the rate of ionization is 1.6x10 and 2.6x10 (H atom sec) ^ 

7 

for cases (1) and (2) , and the rate of Li production (using the cross sections 

"26 “X 

of Figure 1) is 6x10 (H atom sec) for both these cases. Thus, our local 

-15 7 

conditions do not conflict with an ionization rate of .^J-O ; the Li production 

could be reduced by a lower alpha particle abundance in the low energy cosmic 

rays and a steeper energetic particle spectrum. Furthermore, processing of 

7 

interstellar matter through stars could destroy some of the Li, and the 

infall of extragalactic gas could introduce further uncertainties in the 

problem of the origin of this isotope (see Reeves 1978) . 

The calculated 4.44 MeV narrow line intensity for cases (3) and (4) 

in Table 5 are larger by about an order of magnitude than that calculated by 

( 

Meneguzzi and Reeves (1975), mainly due to the enhanced abundance of elements 

\ 

heavier than He along the line of sight. 

Bussard et al. (1978) have related the production of nuclear gamma rays 

to x-ray iron line production, q. resulting from charge exchange of 

0.0 

low energy cosmic rays with interstellar matter. For the abimdances of case 
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(2), and s=4 and = 20 MeV/nucleon, they find that (narrow) /q^ g 

is about 4.5x10 This ratio also applies to case (4), since g depends 

only weakly on the ambient medltm abundances (Bussard et al. 1978 and 

R,' W. Bussard private communication 1978) and both qg g and q^ scale 

linearly with the cosmic fay abundances. Thus from Table 5, the 6.8 keV 

2 2 

flux for case (4) should be ^ 0.1 photons/cm sec rad, or .vX).02 photons/cm sec 
from a longitude interval of IS® around the galactic center. This is con- 
sistent with the upper limit of 0.024 given by Bussard et al. (1978) for 
that region. 

The nuclear line emission discussed above is superimposed on the possible 

bremsstrahlung continuum that is apparently required to account for the 

\ 

difference between the observed high energy gamma rays and the -rf** decay 

spectrum at energies between about 50 and 300 MeV (see Figure 29) . It should 

be noted, however, that unresolved discrete sources (Higdon and Lingenfelter 

1976, Bignami, Caraveo and Maraschi 1978) may also account for some of the 

observed high energy gamma ray emission. 

The bremsstrahlung fltix shown by the dashed line in Figure 29 has a 

numerical value given by 5x10 photons/ (cm^ sec' rad MeV). If the 

electrons that produce this flux have the same spatial distribution as 

q o(cu)/q o(0)» then the local bremsstrahlung eraissivity, q , .and the corres- 
TT TT , n 

ponding local electron intensity, are 

q„(E ) =-10 ^^E photons/ (H atom sec MeV), (27) 

B Y Y 

and 

j^(E) =-1.4 E ^ electrons/ (cm^ sec sr MeV); (28) 

equation (28) is obtained from (27) by using the bremsstrahlung formula 
given by Ramaty and Westergaard (1976)., This local electron intensity, j^, 
joins smoothly at- about 200 MeV with the electron intensity . deduced by 
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Goldstein, Ramaty and Fisk (1970) from the galactic nonthermal radio 
emission. However, a power law extrapolation of to energies below 

about 50 MeV is not required by the gamma ray data; nor is it required by 
the galactic nonthermal radio emission. 

The resulting sums of the nuclear and bremsstrahlung fluxes are shown 

in Figures 30 and 31 for the abundance combinations of cases (3) and (4), 

respectively. Many narrow and very narrow lines are visible, especially 

in Figure 30 where the ambient meditom and low energy cosmic rays> have the 

same composition. The intensities of the strongest lines are given in 

Table 5; case (3) corresponds to Figure (30) and case (4) to Figure 31. 

The width of the 4.44 MeV line is about 110 keV and this line has.no very 

narrow component. The very narrow component of the 6.129 MeV line conta in s 

about 20% of the line photons for our assumed conditions (grains with mean 

radius = 5x10 ^cm containing half the interstellar 0). But the relative 

importance of the very narrow component depends strongly on the fraction of 

0 locked in grains, on the grain sizes, and on the spectrum of the energetic 

.particles, as discussed in Illb. The width of very narrow grain lines is 

determined .principally by differential galactic rotation. In Figures 30 and 

31 all of the photons of the very narrow 6.129 MeV line are in a 5 keV bin, 

and about 80% of the narrow plus very narrow 0.847 MeV line photons are in a 

2 keV bin. Other potentially detectable lines seen in these figures are at 

1.238 MeV from ^^Fe, at 1.369 MeV from ^^Mg, at 1.634 MeV from ^^Ne, at 
28 2 

1.779 MeV from Si, at 1.809 MeV from Mg and at 2.313 MeV and 5.105 MeV 
from N. 

-5 -4 

The 0.511 MeV line fluxes in Figure 30 and 31 are 6x10 and 1.7x10 
2 

photons/cm sec rad, as given also in Table 5. In these figures, all the 
0.511 MeV. photons were' placed in a 2 keV energy bin, as indicated. Such a 
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narrow line width woiild result from positron annihilation in a partially 
ionized gas (Bussard et al. (1979). For annihilation in cold clouds, 
however, these authors find that the FWHM of the 0,511 MeV line is about 
5 keV. 

Although many- of the very narrow gamma ray lines shown in Figures 

30 and 31 may be observable only with high resolution detectors, the 
strongest lines and line con5>lexes should be observable with detectors of 
much lower energy resolution. This can be seen from the nuclear spectra 
shown in Figure 29, which are the same as those shown in Figures 30 and 

31 except that they are averaged over broader energy bins as indicated in 
the figure. Even, with such low energy resolution considerable structure is 
still quite evident. Prominent lines at 0.5 and 4.4 MeV and a broad feature 
between about 1 and 2 MeV can be seen for case (4). For case (3), line 
features at approximately 0.85, 1.3, 1.7, 2.3 and 6.1 MeV can also be seen. 

The nuclear lines resulting from low energy cosmic ray interactions 

can be compared with other diffuse interstellar gamma ray lines which are 

expected from decay of radioisotopes synthesized in supernova and nova 

explosions and from annihilation of positrons produced both by decay of 

such radioisotopes and by other processes. 

The most important of such radioisotopes thought to be synthesized in 

supemovae (Clayton, Colgate and Fishman 1969; Clayton 1973; Ramaty and 

* 2 ^ ' "* * * 

Lingenfelter 1977; Arnett 1977) are A^, which produces a line at 1.'809 MeV, 
60 

and Fe with lines at 0.059, 1.173 and 1.332 MeV. Following Lingenfelter 
and Ramaty (1978), we take the production rate of these isotopes at the 
present epoch to be about 10% of the average rates required to produce over 
the age of the galaxy 0.1% and 1% of the galactic ^^Mg and ^*^Ni abundances, 
respectively. The flux in each of the above four lines is then about 

-5 2 

7x10 photons/ (cm sec rad), comparable to our highest predicted 4.44 MeV flux 
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from low energy cosmic ray interactions. The widths of the A<t and '^Fe 
decay lines are about 0.3% of the line energy, resulting primarily from 
differential galactic rotation. The estimated fluxes per unit energy 
interval in -these lines are therefore quite strong; the 1.809 MeV line 
from radioactivity produced in supemovae could be stronger by about an 
order of magnitude than the line at this energy resulting from cosmic ray 
interactions (Figure 31). 

Accretion from a large, cool companion onto a white dwarf leading to 

thermonuclear runaway in- the CNO cycle has been suggested as the cause of 

novae (e.g. Starrfield et al. 1972). If so, novae could also be sources 

of gamma-ray line emission from the decay of synthesized radionuclei (Clayton 

2? 

and Hoyle 1974). The easiest line to detect should be at 1.275 MeV from Na 

decay with a mean life 3.8 years. Since the estimated nova rate is ^0/year, 

22 

the time between nova explosions in the galaxy is much shorter than the Na 

mean life. Thus, a detector with a broad field of view observing the galactic 

plane in the direction of the galactic center should observe diffuse emission 

22 48 

at 1.275 MeV. For a Na yield per nova of 10 nuclei (Truran, Starrfield 
and Sparks 1978) and the above nova rate, the intensity of the 1.275 MeV 
line should be comparable to the intensities of the and ^^Fe lines 


given below. 

rm- •, 26 , 22 

The decays of and Na are accompanied by positron emission 


with 0.85 and 0.9 positrons per 1.809 and 1.275 MeV photons, respectively. 


Taking, as above, 0.6 photons per positron in the 0.511 MeV gamma ray line 

(because of positroniiaa formation), we find that these decays should produce 

-5 2 

a 0.511 MeV line with intensity of about 7x10 photons/cm sec rad. This is 


comparable to the 0.511 MeV line intensity resulting from low energy cosmic 
rays for cases (2) and (3) , and is about half the value produced by cosmic 
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rays in case (4) (see Table 5) . Positrons in interstellar space could also 

originate from pulsars (Sturrock 1971), and from shorter lived radioisotopes 
44 

such as Ti (Clayton 1973). 

We now cou^are these calculations with the available observations. 

Haymes et al. (1975) have observed gamma ray lines from the galactic 

-4 

center at 0.53 H<).01 MeV and 4.6 +0.1 MeV with fluxes (8.0 +2.3)xl0 and 
-4 2 ■ 

(9.5 +2.7)xl0 photons/cm sec,’ respectively, and a feature between 1.2 to 

-3 2 

2 MeV with flux (2.6, ^.6)xl0 photons/cm sec. The detector had an opening 
angle of 13° and was centered on the galactic center. The measured width 
of the^4.6MeV line was about 0.7 MeV, larger than- the detector resolution. 

Line emission above about 4 MeV should be due to energetic particle 
reactions; other mechanisms, such as decays of long lived radioactive nuclei 
from processes of nculeosynthesis are not expected to produce lines at these 
high energies . For low energy cosmic ray interactions in the interstellar 
medium, the calcialated diffuse fluxes (Figure 29) in a 0.7 MeV energy bin 
around 4.44 MeV and in a longitude interval of 13° are 2.6x10 ^ and 5.8x10 
photons/cm sec, for cases (3) and (4), respectively. 'This assumes that the 
gamma ray intensity is essentially independent of longitude within 1 radian 
of the galactic center, based on the rough constancy of the high energy 
gamma ray intensity in this longitude interval. 

These calculated fluxes are smaller than the observed flux at 
4.6 +0.1 MeV by more than an order of magnitude. Similarly, the calculated 
fluxes in the broad feature between 1.2 to 2 MeV and in- the. 0.51 MeV line 
are also much lower than those observed. But the line ratios, nonetheless, 
are consistent, within errors, with the ratios expected from energetic 
particle reactions. 

To increase the fluxes calculated for low energy cosmic ray interactions 
to values that would match the observations would violate the constraints 
discussed earlier. Thus, it is unlikely that the lines observed by Haymes 
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et al. (1975) are produced by such cosmic rays in the interstellar medium. 

Additional information on the nature of the source of the lines 
observed by Haymes et al. (1975) can be obtained by comparing these obser- 
vations with the data of Gilman- et al . (1978) shown in Figure 29 . From 

-2 

measurements on board Apollo 16, they have set an upper limit of 1.2x10 
2 

photons/cm sec on the gamma ray emission from the galactic disk in the 
longitude interval -50“ to' 22°. This upper limit can be compared with the 
data of Haymes et al. (1975) as follows. 

The measured flux of 9.5x10 photons/cm sec at rJy.6 MeV in the energy inter- 
val ,jQi7 MeV and longitude interval of 13“ (Haymes -.et al. 1975) implies a 

-3 —3 

total nuclear gamma ray flux from 0.6 to 5.2 MeV of 4.0x10 and 9.0x10 
2 

photons/cm sec in the same longitude interval, for cases (3)' and (4) respec- 
* 

tively. These values are lower than the upper limit of Gilman et al. (1978) 

and hence consistent with it as long as the .longitude distribution of the 

nuclear gamma rays is sufficiently strongly peaked toward the galactic center. 

However, if the nuclear gamma rays have the same Fongitude distribution as 

the high energy gamma rays (Fichtel et al. 1975), the fluxes of nuclear 

radiation from 0.6 to 5.2 MeV in the longitude interval -50° to 22“ should 
-2 —2 2 

be about 2x10 and 4.5x10 photons/cm sec for cases (3) and (4), and these 

' —2 2 
are larger than the upper limit of 1.2x10 photons /cm sec (Gilman et al. 

1978). We conclude that the nuclear radiation reported by Haymes et al. 

(1975) is mostly confined to longitudes not exceeding the interval of about 

. r 

13° around the galactic center. The possibility of producing nuclear gamma 
rays by a massive black hole at the galactic nucleus was suggested by 
Lingenfelter, Higdon and Ramaty (1978). 

Leventhal et al. (19 78a, b) have recently made high resolution spectral 
measurements of positron annihilation radiation from the galactic center. 

-3 

The observed line at 0.511 MeV line has an intensity of (1.21 +0.22)xl0 
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photons/cm sec), and a FWHM of less than 3.2 keV. There is also evidence, 

for the 3-photon continuum from triplet positronium annihilation. The 

detector used had an opening angle of 15° centered on the galactic center. 

The observed line intensity is consistent with the lower resolution 

measurements of Haymes et al. (1978) of the line at ~0.53 MeV. 

By the same argument as presented above for the '4.4 MeV line, the 

0.511 MeV line intensity measured by Leventhal et al. (1978a, b) is also 

much greater than that expected for diffuse emission from the interstellar 

medium. The intensity of the 0.511 MeV line from low energy cosmic ray 

-5 -4 '2 

Interactions is, from Table 5, 6x10 and 1.7x10 photons /cm sec rad, for 

cases (3) and (4), respectively. Sinc'e the intensity of this line from 

26 22 '^5 2 

A.fand Na decays is ..^7x10 photons /cm'^ sec rad -5 as discussed above, the 

-4 

total 0.511 MeV line intensity from ’diffuse emission is about 1.3x10 and 
“4 2 

2.4x10 photons/cm sec rad for these two cases'. If we assume that the 

longitude dependence of this emission is the same as that of the high 

> 

.energy gamma rays, as we did for the 4.4 MeV line, then we "dhtain 0.511 MeV 

-5 —5 2 

line fluxes of 3.4x10 and 6.3x10 photons/cm sec in the longitude interval 
of 15° around the galactic center, for the above two cases, respectively; 
these values are much smaller than the observed' flux. 

The origin of the positrons responsible for the observed 0.511 MeV 
line is not known, and all the mechanisms discussed alJove (see also 
discussion by Ramaty 1978) coiild be possible candidates. 

If, however, the positrons are produced predominantly in energetic particle 

reactions, then the 0.511 MeV line should be accompanied by strong nuclear 

12 

lines such as the 4.44 MeV line from G deexcitation. The observations of 
Haymes et al, (1975), discussed above, would support this possibility as we 
shall: show, but these -measurements need to be confirmed by independent 


observations. 
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The relative. intensities of the measured positron annihilation and 

12 

C deexcitation lines are consistent with energetic particle reactions. 

For s=4 and = 20 MeV /nucleon, the ratio of positron production to gamma 
ray production in the 0.7 MeV interval around 4.4 MeV is 0.82 and 1.1, for 
cases (3) and- (4), respectively; for 0.6 gamma rays in the 0.511 MeV line 
per positron, the ratio of the 0.511 MeV line to these gamma rays is 0.5 
and 0.65 for these two cases. The observed ratio is 1.27 ^.42, if we use 
the Leventhal et al. (1978a, b) flux for the 0.511 MeV line, and 0.84 ^.34, 

,if we use the flux given by Haymes et al. (1975) for the line at ~0.53 MeV. 

(The line at 6 MeV was measured only by Haymes et al. 1975). The latter 
ratio is consistent, within errors, with the calculations, while the former 
would allow some additional positrons from sources other than energetic 
particle reactions. But in any event, the positive detection of a 4.4 MeV 
line with intensity comparable to that observed by Haymes et al. (1975) would 
imply that a major fraction of the positrons are produced in energetic 
particle reactions. 

As discussed above, the source of the nuclear radiation' responsible 
for the lines observed by Haymes et al. (1975) should have a small angular 
size «13°) since otherwise there is a conflict between the observed line 
intensities and 'the upper limits of Gilman et al. (1978). A similar argu- 
ment can be made for the source of the 0.511 MeV line if the positrons are 
produced by energetic particle reactions. Moreover on the basis of the 
very narrow observed width of the 0.511 MeV line, Bussard et al. (1979) find 
that a large fraction of the positrons should annihilate in an ionized 
medium,' and a possible implication of this result is that the positrons 
originate in a point-like source at the galactic center and annihilate 
in the HII regions close to the galactic nucleus. The size of the" 0.511 MeV 
line source would in this case be less than about 2“. Johnson and Haymes (1973) 
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have observed the galactic center during a lunar occultation and found a 
2.3 standard deviation decrease of the 0.44 to 0,54 HeV flux during the 
occultation. This result also suggests a small (^.5**) angular size for 
the source of .the 0.511 MeV line. 

If the bulk of 0.511 MeV line is from nuclear reactions of energetic 
particles, we can calculate from Table 4 the rate of energy depositions 
associated with the production of the positrons. This amounts to approxi- 
mately 1 erg per 0.511 MeV photon for cases (3) and .(4), so that. the 

43 

observed flux requires' the deposition of .^<10 erg/sec. For a source size 

less than about 2 ° 'this energy deposition rate is larger than the infrared 

luminosity of the galactic center (Hofmann, Frederick, and Emery 1971) by 

about ah order of maghitud'e. Some of the excess energy could be in the form 

of mass motions (Oort 1977) or lost beyond the event horizon of a massive 

black hole (Lingenfelter et al. 1978). Less energy deposition, however, is 

required if the positrons are produced from radioisotopes S3mthesized in 
44 

supernovae (e,g. Ti) or by pulsars. The key to the question of the origin 
of the positrons is the detection of other gamma ray lines Jthat are character- 
istic of the prevailing positron production. mechanism. 
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V. SUMMARY . 

We have studied the production of gamma ray lines from the interaction 

V 

of energetic 'particles with the abundant constituents of cosmic matter. 

We consider the reactions induced by energetic protons and alpha particles 
in ambient nuclei (He, C, N, 0, Ne, Mg, AJi, Si, S, Ca and Fe) , and the 
inverse reactions in which energetic nuclei interact with' ambient H and 
He. We have evaluated the line production cross sections by analyzing 
a large body of laboratory nuclear data which we discuss in detail. We 
consider prompt gamma rays produced by direct excitation of nuclear levels 
and by spallation reactions which leave. the secondary nucleus in an excited 
state, ,as well as delayed emission resulting from long-lived radioactive 
nuclei that are also produced in the energetic particle reactions. A list 
of all the lines that we consider is given in Table 1; in addition, our 
calculations also take into account unresolved nuclear radiation from targets 
heavier than 0. The nuclear data is discussed systematically in Section 
III. The lines of He are treated in Ilia, those of C, N and 0 in Illb, 
while those of Ne, Mg, A£, Si, S, Ca and Fe are presented in IIIc. 

We detemine the shapes of the gamma ray lines by taking into account 
the kinematics of the reactions, data on the angular distribution of both 
the secondary particles and the gamma rays, and the, lifetimes of the excited 
levels or their radioactive progenitors. If these lifetimes are sufficiently 
long, very narrow lines can be produced in solid targets such as interstellar 
grains. The evaluation of line shapes is treated in Ilia, Illb and IIIc, 

We have constructed a computer program that evaluates, both digitally 
and graphically, nuclear gamma ray spectra from about 0,1 to 8 MeV, The 
bulk of the nuclear radiation from energetic particle reactions is in this 

I 

energy range, although some nuclear lines of astrophysical interest can 
also be found at higher energies (Crannell, Crannell and Ramaty 1979), 
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In. this computer program it is possible to vary the abundances of both 

the ambient medium and the energetic "particles, the energy spectrum and 

angular distribution of these particles, and the amount of matter in grains, 

as well as the characteristic grain radii. Various numerical results based 

on this program are presented in Illd. For comparison purposes, in this 

subsection we also present results on positron production in energetic 

particle reactions based on updated values of our earlier cross section 
* 

compilations (Ramaty et al, , 1975) . 

As an application of the theory developed in this paper, we have 

considered gamma "ray line production in the interstellar medium. The theory, 

however, can also be applied to other astrophysical sites such as solar 

flares, compact objects and nuclei of galaxies. We have evaluated the 

nuclear gamma ray emission from the galactic disk by normalizing the 

emissivity of the nuclear lines to that of high energy gamma rays, by 

using the observed intensities of these gamma rays from the galaxy, by assuming 

that the high and low energy cosmic rays have the same galactic distribution, 

and by .taking into account the possible enhancement of heavy elements toward 

3 

the galactic center. For a local low energy cosmic ray density of 1 eV/cm 

and heavy eluent abundance gradients consistent with x-ray absorption 

-5 -5 -5 

observations, we find intensities of about 7x10 , 3x10 and 10 

photons/cm sec rad for the strongest nuclear lines at 4.44, 6.129 and 
0.847 MeV. The intensity of the 0.511 MeV line from positron annihilation 
for the same assumptions is 10 ^ photons/cm^ sec rad. These lines should 
be detectable by planned high sensitivity observations. Such detection 
could give a wealth of Information on the properties of the low energy 
cosmic rays and the composition and spatial distribution of interstellar 
gas and dust. 
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We have compared the results of our calculations with line observations 
from the direction of the galactic center. The reported fluxes at 0,511 
MeV (Leventhal, MacCallum and Stang 1978 a,b) and at ''j 4.4 MeV (Haymes et 
al. T975) from a region of angular size < 15° around the galactic center 
are larger by more than an order of magnitude than those expected from 
diffuse low energy cosmic ray interactions. By using the theory developed 
in this paper, we find that the sources of these lines should be much 
more sharply peaked around galactic longitude & ~ 0 than those of the 
high energy gamma rays; otherwise, the broad band nuclear radiation that 
accompanies the spectral lines would exceed the upper limits (Gilman et 
al. 1978) on the total 0.6 to 5,2 MeV radiation from the galactic disk 
in the longitude interval -50° to> 22°, It is conceivable that the positrons 
and nuclear lines could be produced by one or more point sources close 
to the galactic center. 

Note added in proof : After the completion of this manuscript, preliminary 
results from the A-4 experiment on HEAO-1 became available (J. L. Matteson, 
private' communication *^1978) . These results do not confirm the Haymes et al. 
0-975) observations of the 4.44 MeV line from the galactic center. The 
upper limit set by the HEAD. observations on the total 4 to 7 MeV radiation 
implies that not more than about 30% of the positrons required to account 
for the observed 0.511 MeV line (Leventhal et^al. 1978a, b) could result 
from energetic particle reactions if the nuclear reaction rate does not vary 
on a time scale shorter than the positron annihilation time. 
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Figure Captions 

1. Cross sections for alpha particle reactions with He. Solid 

curve - ^Be production in its ground state plus excited state. 

7 • 

Dashed curve. - Li production in its excited' state. 

• . 12 

2. Gross sections for 4.44 MeV photon production from C. Solid 

and dashed-dotted curves - direct excitation of the 4.439 MeV 
level in (p>p^) and (of, c/) reactions. Dashed and- dashed-crossed 
curves - total production of 4.44 /MeV photons by protons and 

alpha particles, respectively. 

( 

3. Gamma ray production cross sections in spallation reactions of 

X2 X6 

-protons and alpha particles on G and 0. 

4. Cross sections for 2.313 MeV photon production from and ^^0. 

Solid curve - 2. ,313 MeV photon production cross sections hy direct 

14 

proton excitation of levels in N; dashed curve total 2.313 MeV 

14 14 

photon production by protons on N including 0 decay; dashed- 

crossed curve - 2.313 MeV photon production by direct alpha 

14 

particle excitation of levels- in -N; dashed-dotted curve - total 

X6 

2.313 MeV production from proton bombardment of 0. 

5. Gross sections of a-/ 6; 2 MeV photon production from ^^0. Solid 


and dashed-dotted curves - 6.129 MeV photon production by direct 
proton and alpha particle excitation of levels in ^^0. Dashed 


and dashed crossed curves - 6.1 to 6.3 MeV photon production from 

X6 

proton and alpha particle bombardment of 0. 
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6. Cross sections for 6,917, 7.117, 7.299 and 2.741 MeV photon produc- 

”16 “ 

tion by protons and alpha particles on 0. The 6.917, 7.117 and 

16 

2.741 MeV lines are from deexcitation of levels in 0, and the 

15 ■ 

7.29-9 MeV line is from a level in the spallation product N. 

20 

7. Gamma ray production from Ne. Solid and dashed— dotted curves — 
cross section for the direct excitation of the 1.634 MeV level 

in Ne by (p,p^) and (o:,Qf^) reactions. Dashed and dashed-crossed 

curves- total 1,634 MeV photon production from proton and alpha 
20 

particles on Ne; these cross sections include the contributions 

20 

of cascades from higher lying levels in Ne. Crossed and crossed- 

dotted curves - total photon production from proton and alpha 

20 

particle bombardment of Ne (Ey> 0,7 MeV). 

24 

8. Gamma ray production from .Mg. The various curves represent 

quantities similar to those given by the cruves of Figure 7. 

■> 

2 6 ' 22 

9. Cross sections for the production of A-L and Na by protons on 

isotopes of Mg and Si. The dashed curve for proton bombardment 

2^ 23 26 

of Mg includes the contributions of Mg, Mg and Mg, using 

the isotopic ratios of Cameron (1973). 

28 

10. Gamma ray production from Si. The various curves represent 

quantities similar to those given by the curves of Figure 7 . 

56 

11. Gamma ray production from Fe. The three curves represent 


quantities similar to those given by the curves of Figure 7. 


Gamma ray production by alpha particle bombardment of Fe appears 
to be small at low energies because of the high Coulomb threshold. 
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12 . 


13. 


/V 

u. 


15. 


Cross sections for ^^Co, ^^Co, '^^Mn and ^^Mn production by- 
protons on ^^Fe. 

' 12 16 

The angular distribution of recoil C and 0 in the center of 

mass frame for the reactions ^^C(p,p' )^^C* and ^^0 (pjp'^ ) ^^0 . 

The profile of the 4.44 MeV gamma ray line from the reaction 
12 12 

C(PjP^) C*. The bombarding protons are monoenergetic and 
confined to a beam, and the gamma rays are observed at 90® to the 
beam. The top curves are calculated for an uneven population of the 
spin states of the level. The lower curve is 

calculated for isotropic gamma ray emission. 

The profile of the 4.44 MeV line gamma ray line from the reaction 
^^C(p,p' )^^C* for an isotropic distribution of bombarding 
protons with an energy spectrum proportional to E“^, The top curve 


is for the case of uneven spin state populations, and the lower 
curve is for isotropic gamma ray emission. 
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16. The mean stopping ranges, <x>, the energy loss rates, dE^/dx, and 

X 6 *2J\ 5 6 

the range fluctuations, for low energy 0, Mg and °Fe 

stopping in water. 

- 1-7-i The energy distribution in the laboratory frame of unresolved 
gamma rays from proton bombardment of heavy nuclei. 

18. Monte Carlo simulated gamma ray spectrum for energetic particles 
and ambient medium having solar compositions; s and E^, are the 
spectral parameters of the energetic particles, and is the 
characteristic radius of the interstellar grain distribution. “ 

19. Monte Carlo simulated gamma ray spectrum for energetic particles 
and ambient medium having solar compositions . 

20. Monte Carlo simulated gamma ray spectrum for energetic particles 
and ambient medium having solar compositions-. 

21. Monte Carlo simulated gamma ray spectrum for energetic particles 
having cosmic ray source composition and ambient medium with solar 
composition. 

22. Monte Carlo • simulated gamma ray spectrum for energetic particles 
having cosmic ray sourc.e composition and ambient medium with solar 
composition. 

23. The gamma ray spectrum of Figure 21 averaged over energy intervals 
of 100 keV. Dashed curve — .photons produced by only nuclei 
heavier than He; solid curve - photons produced by all the cosmic 
rays . 

The gamma ray spectrum of Figure 22 averaged over energy intervals 
of 100 keV. Dashed curve - photons produced by only nuclei heavxer 


24 . 
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than He; solid curve - photons produced hy all the cosmic rays . 

25. Gamma ray emissivity per unit energy density in the energetic 

particles. The emissivities are for the narrow 4.44 MeV line, 

the narrow plus very narrow 6.129 and 0.847 MeV lines, and the 

total broad, narrow and very narrow emissions. Both the energetic 
particles and ambient ;9ediui]; have solar compositions ? The parameters 
s and Ec are- defined in equation (22). 

26. Emissivities of the narrow 4.44 MeV and narrow plus very narrow 
6.129 MeV lines per unit energy deposition rate of the energetic 
particles in a neutral medium. Both the energetic particles and 
ambient medium have solar compositions . 

27. Emissivities of the narrow 4.44 MeV and narrow plus very narrow 
6.129 MeV lines per unit energy deposition rate of the energetic 
particles in a ionized medium. Both the energetic particles and 
ambient medium have solar compositions. 

28. Ratios of narrow plus very narrow 6.129 MeV photon production, and 
positron production, to the narrow 4.44 MeV production. The 
ambient medium has solar abundance, and cy/p = 0.1 in the energetic 
particles . 

29. Gamma ray emission from the interstellar medium from the direction 
of the galactic center. The tt'’ decay spectrum has been calculated 
by Stacker (1976) and is normalized here to the data of Paul et al. 
(1978) . The nuclear gamma ray spectra are calculated for a low 

'energy cosmic ray component with s = 4, E^, = 20 MeV/nucleon, local 

3 ' 

energy density 1 eV/cm , and the same spatial gradient as that 
implied for the high energy cosmic rays from the high energy gamma 
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30.. 


31. 


ray data. The compositions of the ambient medium and low energy 
cosmic rays are for cases (3) and (4) (See text). The dashed line 
is a siB 5 )le power law that could arise from electron hremsstrahlung. 
-High .resolution plot of the. sum of the nuclear gamma ray and bremsstrah- 
lung fluxes from the direction of the galactic center. The parameters 
are the same as for Figure 29 case (3) . 

High resolution plot of the stm of the nuclear gamma ray and bremsstrah- 
lung fluxes from the direction of the galactic center. The parameters 
are the same as for Figure 29 case (4) . 
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TABLE 2 


^ - a 


NUCLEUS 

ARTTWnANCF. 

DIRECT DEEXCITATION 
LINES (MeV) 

SECONDARY NUCLEI AND THEIR 
GAMMA RAY LINES (MeV) 


1 

— 


^He 

10“^ . 

— 

®Li(3.561) ,^Li(0.478), 
^Be(0.431;0.478) 


4.2x10“^ 

4.438,15.11 

^®B(0.717 ;1.0p) ,^°C(0.717) , 
^^B(2.124;4.443;6.741;6.791) , 
^^C(1.995;6.337;6.478) 


4.7x10"^ 

3.684,3.853 

— 


8.7x10"^ 

1.632,2.313,5.105 

^^C(4.438;15.11) ,^S(2.313) 


6.9x10“'* 

2.741,6.129,6.917,7.117 

^°B(0.717),^°C(0.717), 

^^C (4. 438; 15. 11) , ^®C(3.684; 
3.854), ^^N(1.632; 2.313; 5.105), 
^^0(2.313) , ^®N(5.270;5.298;6.322; 
•7.299), ^^0(5.180;5.241;6.176) 

20 

Ne - 

3.JxlO“^ 

1.634^2.613,3.334 

^®0(6.129) ,^^F(0.110;0.197-)., 

19 20 

Ne(0.238;0.275) , ''Na(i:634) 


4.1x10“^ 

1.275 

22 

Na(1.275) 


3.2x10“^ 

1.369,2.754 

20 ’ 22 
Ned. 634, 2. 613) , Na(0.891;1.275), 

^®Na(0.440;1.636;2.640) , 

^®Mg (0.451; 1.600) 

25„ 

Mg 

4.0x10“® 

— 

^^Na(1.275) 


4.4x10“® 

1.809 

^®A1 (1.809) 

2^A1 

3.3x10“® 

1.014 

^®A1(1.809) 


4.5x10“® 

1.779,5.099,6.878 

^®Ne (1.634; 2. 613) ,^^Na(1.275) , 
^'*Mg(l.‘369) ,^®A1(1.809) , 
^^Al(0.844;1.0l4);^''si(0.780;0.957 


1 . 6xl0“® 

2.230 

^®Si(1.779) ,®^P(1.266;2.029; 
2.234), (1.249; 2. 034; 2. 232) 

,'°Ca 

2.2x10“® 

3.736 

— 


3 . 2x10“® 

0.847,1.238,1.772, 

®^Cr(0.744;1.334;1.434) 



1.811,2.094,2.113 

®^Mn ^0.744; 0.936; 1.434) 
®^Mn(0.835), ®^Fe(1.408) 
®®Fe(0.092;0.412;0.477;0.931; 
1.370;1.408), ®®Co(0.477;0.931; 
1. 408), ®®Co(0.158;0. 812; 0.847; 


. -1.038;1.238;1.772;2.599) 



TABLE 3 


ABUNDANCES 


Element Solar System 


Local 

Cosmic Ray Source 


5kpc 

Ambient 


5kpc 

Cosmic Ray Source 


H 

1 

. 1- ■ 

^ - 1 

1 

He 

10-1 

0.1 

0.1 

0.1 

C 

4.2x10“'^ 

3. 8x10“ 1 

2.1xl0"l 

1.9x10-^ 

N 

8.7x10”^ 

5 . 6x10 ^ 

1.3xl0“l 

8.4xl0“l 

0 

■ 6.9x10“^ 

4.3xl0"l 

3. 5xld"l . 

2 .-2x10-^ 

Ne 

3^7x10“^ 

S.SxlO"*^ 

1.9x10"'^ 

2.9xl0“l 

Mg 

4.0x10"^ 

9 .2x10-^ 

2.0x10"'^ 

4.6x10"! 

A-t 

3.3x10"^ 

5.9x10"^ 

1.7x10“^ 

3.0x10“^ 

Si 

4 . 5x10 ^ 

8.1x10“^ 

2.3x10“^ 

4.1x10"! 

S 

1.6x10“^ 

1.2x10“'^ 

8.0x10"^ 

6.0x10 ^ 

Ca 

2.2x10“^ 

9.0x10“^ 

1.1x10“^ 

4.5x10“'^ 

Fe 

3.2x10“^ 

8.5x10-^ 

1. 6x10"^ 

4.3x10"! 



TABLE 4 


4.44 MeV Gantma Ray Line and Relative Positron Emissivities 
Per Unit Energy Content and Deposition Rate of tlie 
Energetic Particles for s=4 and E(.=30 MeV/nucleon 


. <q(4.44)/Vfc> <q(4.44)/Vt> 

[Photons/ (sec Hatom eVcm. ;] [Photons/erg] 



Narrow 

Broad 

Narrow 

Broad 

CASE 

Line 

Line' 

Line 

Line 

(1) 

2.6x10"^^ 

2.6x10"^^ 

0.29 

O'. 29 

(2) 

2.3x10“^^ 

-25 

1..7xl0 ^ 

0.15 

1.1 

(3) 

1.5x10"^^ 

-25 

1.5x10 

1.3 

1.3 

(4) 

8.4x10“^^ 

6.4x10 

0.30 

2.3 


ifc 

q (4.44) + q/ , (4.44 ) 

^narrow ^broad 

(positrons/ photon) 

0.63 

0.83 

0.63 

0.77 



Table 5 


Gamma Ray Line Fluxes from the Direction of the 
. Galactic Center for s=4 and E^.=20 MeV/nucleon 
(in photons /cm2 sec rad) 


.CASE - - 

■ '^4‘/44 (narrow)'- 

?6 ..i.29'(narrow) 
ig’ 129 (very narrow) 

'^dV84 7 (narrow) + 

Iq 847^^®^^ narrow) 

^0.511 

1 

1.3x10“^ 

6.4x10“^ 

2.0x10"® 

lxl0“® 

2 

1.2x10“'^ 

5.9x10“® 

1.9x10"® 

5x10“^ 

3 

7.2x10“^ 

3.0x10“^ 

9.2x10"® 

6xl0“® 

4 

4.2x10“^ 

1.8x10“^ 

5.4x10“® 

1.7x10"'^ 
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Fig, 3 










E (MeV/NUCLEON) 











ORIGINAL PAGE IS 
OE POOR QUALiry 



E(MeV/ NUCLEON) 







24., , ,.24^, *1.369 

Mg(p,p) Mg 

— ^V(p.pV|3e/‘'Mg 
24 

XXX Mg(p,xr) 


24,, , .,24„ *1.369 

— •— Mg ( a,a ) Mg 

— X-^"^Mg (a.a’y ,-J^'^Mg 


cx MgCp.xy) *x*x Mg(a,xy) 

▲ ESTIMATED FROM 5 £ ZOBEL et al (1968) 

ZOBEL et al (1968) ’ .r 

. A FOLEY et al (1962b) 


xX>4*< 


x^ 

X 

X 

x 

X 

X 

X ^ 

r\ ""x 


^Xx T • 

tXx„ . X. 


♦(< 

X, X 

X X 


EMeV/ NUCLEON) 






cr{mb) 




28os/_ _i \28^,*i.779 t i i. f / » 

Si(p,p ) SI — . Si(a,a ) Si 

--- re®i‘'’-P'7.779)''si -X-^/si(a.aV,^„)^«Si 
XXX Si(p.xr) -X-X- Si( 0 ,xr) 


28^. *1.779 


FOLEY et al (!962b) 
^ ESTIMATED FROM 




X X 

• X 

X X 

X* ,x* 

x' i 

• X __ 

X X 

• X _ 

\ 

?,X-X. / \ \ 


'•x. ZOBEL et al (1968) 
’‘•x 

^*x *x. 

X X 

^Xy *X»X»x»X»X»X« 

^^xxxxxxxxxxxx: 


E MeV / NUCLEON 







A 


E 


— ^®Fe(p,p') 0.847 

--^®Fe(p,pVo_847) 


xxxxxxxxxxx 


^ ZOBEL etal. (1968) 

^ ZOBEL et al. (1968) 

A CHANG etal. (1974) 

o JASTRZEBSKI etal. 
(1976) 

^ A ESTIMATED 
\ FROM ZOBEL 
V et al. (1968) 


2 





cr (mb) 



Fig. 12 



l/o-(E)] / (do- /da*) da* 



Fig. 13 





PHOTONS/ (MeV SR) 








ENERGY DISTRIBUTION 



Ey (MeV) 


Fig,. 17 













PHOTONS /keV 



Fig. 20 







PHOTONS /keV 



Fig. 22 









q (PHOTONS/SEC) . q (PHOTONS/SEC) 



0 10 ■ I 10 ,10 


Ec (MeV/NUC) Ec (MeV/NUC) 


Fig. 25 


q/W (PHOTONS /ERG) 



Fig. 26 


q/W (PHOTONS/ERG) 





NUCLEON) 


28 


PEUfilNAL PAGE lu 
E.OOE QUALI'TY 






Fig. 30 




FLUX (PHOTONS CM SEC RAO MeV ) 



Fig. 31 



BIBLIOGRAPHIC DATA SHEET 


1. Report No. 2. Government Accession No. 

TM-79654 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Nuclear Gamma Rays From Energetic 
Particle Interactions 

5. Report Date 

November 1978 

6. Performing Organization Code 
660 

7. Author(s) 

R. Ramaty, Benzion Kozlovsky, R. Lingenfelter 

8. Performing Organization Report No. 

9. Performing Organization Name and Address 

Laboratory for High Energy Astrophysics 
Theoretical Group 

10. Work Unit No. 

1 1 . Contract or Grant No. 

13. Type of Report and Period Covered 

1 2. Sponsoring Agency Name and Address 

14. Sponsoring Agency Code 

15. Supplementary Notes 

To be published in Astrophysical Journal 
Supplement (1979) 

16. Abstract 


Gamma ray line emission from nuclear deexcitation following ener- 
getic particle reactions is evaluated. The compiled nuclear data and 
the calculated gamma ray spectra and intensities c.an be used for the 
study of astrophysical sites which contain large fluxes of energetic 
protons and nuclei. A detailed evaluation of gamma ray line production 
in the interstellar medium is made in the present paper. 


17. Key Words (Selected by Author(s)) 18. Distribution Statement 


Gamma ray lines, nuclear reactions, 
interstellar medium, galactic center 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* ' 

UN 

UN 




For sale by the National Technical Information Service, Springfield, Virginia 22151. 


eSFC 25-44 (10/77) 





















